
Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


2000-09 

Parabolic equation modeling of bottom 

interface and volume reverberation in shallow water 

Li, Lit Siew. 

Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/9409 


Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


http ://w w w. nps.edu/lEbrary 


Calhoun is the Naval Postgraduate Schools public access digital repository for 
research materials and institutional publications created by tire NPS community, 
Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed — and published — scholar^ author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 



NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



THESIS 


PARABOLIC EQUATION MODELING OF BOTTOM 
INTERFACE AND VOLUME REVERBERATION IN 


SHALLOW WATER 


by 


Li, Lit Siew ^ 


September 2000 

Thesis Advisor: 

Kevin B. Smith 

Second Reader: 

James V. Sanders 


Approved for public release; distribution is unlimited. 


duc 


20001222 098 





REPORT DOCUMENTATION PAGE 

Form Approved OMB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instruction, searching existing 
data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate 
or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information 
Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction 
Project (0704-0188) Washington DC 20503. 

1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

September 2000 Master’s Thesis 

4. TITLE AND SUBTITLE 

Parabolic Equation Modeling of Bottom Interface and Volume 
Reverberation in Shallow Water 

5. FUNDING NUMBERS 

6. AUTHOR(S) LI, LIT SIEW 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Naval Postgraduate School 

Monterey CA 93943-5000 

8. PERFORMING 
ORGANIZATION 

REPORT NUMBER 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the 
official policy or position of the Department of Defense or the U.S. Government. 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Approved for public release; distribution is unlimited. 

12b. DISTRIBUTION CODE 


13. ABSTRACT ( maximum 200 words) 

A reverberation model based on the parabolic approximation is developed that includes sediment interface 
and volume perturbations. A multipe forward/single backscatter approximation is made, and the structure of 
the solution is found to depend on the two-way propagation with a scattering strength scaling dependent on the 
local properties of the perturbation. The model is implemented for continuous wave (CW) signals to predict 
mean reverberation pressure levels and for broadband pulse signals to generate complex reverberation 
structures in the time-domain. The spatial correlation and statistical properties of these predicted signals are 
then analyzed in an attempt to extract information on the underlying characteristics of the perturbation. 
Preliminary analysis suggests that reverberation due to the volume perturbations decorrelates more rapidly 
over depth than the reverberation due to interface fluctuations, although the differences appear small. 
Additionally, the statistical character of the reverberation structure due to the interface appears as a relatively 
flat spectrum, while the spectrum of the volume reverberation tends to appear colored. Attempts to correlate 
these characteristics with the structure of the perturbations is ongoing. 


14. SUBJECT TERMS: Shallow-water sound propagation, reverberation, bottom 
interface, bottom volume, spatial perturbations, sound speed perturbations, peak 
correlations, time/temporal series, range series, MMPE, MMPEREV 

15. NUMBER OF 
PAGES 142 

16. PRICE CODE 

17. SECURITY CLASSIFI¬ 
CATION OF REPORT 
Unclassified 

18. SECURITY CLASSIFI¬ 
CATION OF THIS PAGE 
Unclassified 

19. SECURITY CLASSIFICA¬ 
TION OF ABSTRACT 
Unclassified 

20. LIMITATION OF 
ABSTRACT 

UL 


NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 


Prescribed by ANSI Std. 239-18 298-102 


1 

















11 




Approved for public release; distribution is unlimited. 







ABSTRACT 


A reverberation model based on the parabolic approximation is developed that 
includes sediment interface and volume perturbations. A multiple forward/single 
backscatter approximation is made, and the structure of the solution is found to depend 
on the two-way propagation with a scattering strength scaling dependent on the local 
properties of the perturbation. The model is implemented for continuous wave (CW) 
signals to predict mean reverberation pressure levels and for broadband pulse signals to 
generate complex reverberation structures in the time-domain. The spatial correlation and 
statistical properties of these predicted signals are then analyzed in an attempt to extract 
information on the underlying characteristics of the perturbation. Preliminary analysis 
suggests that reverberation due to the volume perturbations decorrelates more rapidly 
over depth than the reverberation due to interface fluctuations, although the differences 
appear small. Additionally, the statistical character of the reverberation structure due to 
the interface appears as a relatively flat spectrum, while the spectrum of the volume 
reverberation tends to appear colored. Attempts to correlate these characteristics with the 
structure of the perturbations is ongoing. 
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I. INTRODUCTION 


In 1993, the Office of Naval Research (ONR) sponsored a reverberation program, 
the Acoustic Reverberation Special Research Program (ARSRP), to study the primary 
causes and nature of acoustic reverberation in the deep ocean. The region studied was 
near the Mid-Atlantic Ridge, an area where bottom characteristics are dominated by 
highly variable topography and hard, basaltic rock structures (Tucholke, et al, 1993a and 
1993b). The dominant feature contributing to long-range reverberation was then the 
water/bottom interface. Furthermore, the general structure of the reverberation signal was 
found to correlate quite well with the forward propagation (Smith, et al., 1996, and 
Makris, et al., 1994), i.e. high reverberation returns were found to occur where energy 
interacted strongly with the bottom. 

Previous analysis with a broadband parabolic equation (PE) model has shown 
good agreement with predicted arrival time structures due to deep-ocean bottom interface 
reverberation (Smith and Cushman, 1997). Application of this approach for both single 
frequency (CW) and broadband pulse propagation calculations showed good agreement 
with measured data. In addition, the broadband reverberation predictions indicated that 
some features of the reverberation signal required the discrimination of multipath 
propagation in order to fully resolve the structure, particularly in CW shadow zones. 

A reverberation experiment is now being organized by ONR in the littoral regions 
of the East and South China Seas called ASIAEX. In contrast to the previous deep ocean 
regions studied, this area has a much smoother and softer water/bottom interface. It is 
unclear at this point whether the dominant mechanism for bottom reverberation will be 
interface roughness or volume inhomogeneities in the sediment. One may expect that the 
measured reverberation will contain statistical information related to the character of the 
scattering mechanism. However, the reverberation signal from the sediment volume may 
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be affected by forward scatter from the interface and vice-versa, thereby complicating the 
distinction of the dominant mechanism. 

The purpose of the present work is to establish a model for predicting the 
influence of propagation on both interface and volume reverberation in a shallow water 
environment. This model is similar to the previous PEREV model (Tappert and Ryan, 
1989), although it formally treats the scattering influence in the context of the PE 
approximation. It is further expanded to treat volume scattering. 

Starting with based on the Monterey-Miami Parabolic Equation (MMPE) 
propagation model developed by Smith and Tappert (1994) and Smith (2000), bottom 
interface and sediment sound speed perturbations were incorporated and broadband and 
CW analysis conducted. The theoretical basis of reverberation due to the bottom interface 
and bottom volume is developed in Chapter II. Specifically, we were able to use a CW 
one-way propagation model to show that the reverberant field is the coupling of two, one¬ 
way propagations. The theoretical basis for reverberation was expanded to time-domain 
analysis for both the interface and volume in order to provide simulation for pulsed 
propagation. In this method, it was shown that the two, one-way propagations from the 
source to the scattering patch and from the receiver to the scattering patch would give the 
reverberation field by convolution of both fields in time. Chapter II also details the 
models used for generating realizations based on statistically controlled random 
fluctuations for the interface and volume, and offer a practical implementation that was 
incorporated into the MMPE model. This new model is known as MMPEREV (see 
Chapter HI). 

Secondly, after incorporating random fluctuations of pertinent parameters into the 
MMPE model, post-processing analysis of these results in MATLAB was performed. 
Chapter IV describes the post-processing and signal processing implemented in 
MATLAB and provides a visual presentation of the effects of reverberation. A single set 
of environmental parameters was used in all the simulations to afford a controlled 
computation of the effects of reverberation, providing a qualitative description of the 
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analysis. Comparisons between perturbed and unperturbed data were conducted and 
detailed in Chapter IV. Post-processing and signal analysis included: 

• time-domain analysis and comparison with unperturbed results; 

• CW analysis and comparison with unperturbed results; 

• temporal and range peak auto- and cross-correlations in depth; 

• spectral analyses of the reverberation and reverberation loss. 

This model will serve as a tool for understanding propagation influences in future 
reverberation experiments. 



THIS PAGE IS INTENTIONALLY LEFT BLANK 
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II. THE REVERBERATION PROBLEM - SHALLOW WATER 
REVERBERATION AND SCATTERING THEORY 


In this chapter, we present a theoretical treatment of the effects of bottom 
reverberation and scattering. We acknowledge that in shallow water, the effects of 
reverberation are more involved due to the thin layer of water and the somewhat invariant 
sound speed profile in the water column. However, perturbations to the interface and 
sediment volume (hereafter known as just the ‘volume’) exist and may be modeled 
statistically. These will be treated in detail. One element lacking at this point is the study 
into density variabilities in the volume. Implementation for density variability will form 
part of the scope of follow-on work, which is discussed in Chapter V. 

The interface and volume perturbations may be treated separately because of the 
nature of their perturbations - the interface primarily being perturbed spatially while the 
volume variabilities involve both the sound speed and density. In order to effectively 
model their total effects, all interface and volume perturbations must be jointly 
incorporated into the model at the same time. This has been integrated into the MMPE 
Model developed by Smith, 2000 (originally known as the UMPE Model developed by 
Smith and Tappert, 1994). 

A. BOTTOM INTERFACE SCATTERING 

We begin with the standard wave equation, also known as the Helmholtz wave 
equation. It is written for acoustic pressure, p, as 

V 2 p + k 2 n 2 p = 0 (2- 1 ) 

with the following definitions for wavenumber and index of refraction, including the 
reference sound speed, 

k 0 =O)lc 0 and n(x) = c 0 / c(x). ( 2 - 2 ) 
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Now, consider a rough bottom where the spatial fluctuations may be defined 
statistically. The boundary between the water column and the bottom is then defined by, 

n 2 =n 2 w (l-H(z-z b )) + n 2 b H(z-z b ), (2.3) 

where n w and n b represent the spatially varying indices of refraction in the water column 

and bottom (volume), respectively, and z b is the bathymetric depth at the interface. The 
definition of the Heaviside step function used at the water/bottom interface is given by 


H(z) = 


' 2 ’ 


0 , 


z >0 

z=0 , and ^-H(z) = S(z). 


(2.4) 


We treat the influence of the roughness of the bottom interface by introducing 
perturbations to the bottom depth of varying length scales. This is akin to the method of 
small perturbations (MSP) found in Brekhovskikh and Lysanov (1990). More 
specifically, we define 

Z b (r) = z b (r) + tj l (r) + T] s (?), (2.5) 

where z b (?) is the deterministic portion of the bathymetry as a function of horizontal 
range, F. The deterministic part of the bathymetry may also be seen as the mean depth on 
a relatively flat portion .of the seabed. Correspondingly, tj t (?) is the long-wavelength 
portion of the randomly generated roughness comprising a spectrum of roughness, while 
7j s (r) is the short-wavelength component, also generated randomly. The short- 

wavelength components are assumed to be the components generating scatter, thereby 
responsible for reverberation, predominantly at the Bragg wavenumber. The short 
wavelength components may be thought of as the easily observeable variability of the 
seafloor, i.e. loose impediments and small mounts. The longer wavelengths are then the 
larger undulations. Hence, it is reasonable to think of a ‘critical’ or ‘limiting’ 
wavenumber separating these two components. This critical or limiting wavenumber is 

assumed to be approximately, k Um ~^-k 0 d*, where 0 C is the critical angle of 
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propagation, and smaller roughness wavenumbers (i.e. longer wavelengths) are assumed 
to dominate the forward scatter properties of the acoustic field. 

Performing a Taylor series expansion of the Heaviside step function we have 
defined in Eq. (2.4), we can then re-write the original Helmholtz equation, Eq. (2.1), as 

V 2 p + k 2 (n 2 +fi)P = 0 (2.6a) 

where 

n 2 =nl(l-H(z-z b ~rii))+n 2 b H(z-z b -ri l ) (2-6b) 

and 

ju(r, z) = (n 2 w - n 2 )r] s 8(z -z b -r)i)- ( Z6c ) 

This shows that the primary part of the propagation is due to the long-wavelength 
components and the perturbation is due to the small wavelength components. The small- 
scale features previously expressed as TJ s (r) are considered to introduce secondary cross¬ 
terms whose contribution will be small and may thus be ignored. From hereon, ignoring 
the short-wavelength components will introduce an approximation to the theoretical 
development and variables will then be “hatted” to denote the use of only long- 
wavelength perturbations, e.g., the approximate refractive index, n. 


We now express the wave equation in cylindrical coordinates as 

Ip + (n 2 + n)p = 0. 


~ld_ 

d ) 

i a 2 a 2 ' 

r dr 

1>J 

l + r 2 a ^ 2 + az 2 _ 


(2.7) 


Now let p = —, then 
Vr 


d 2 u . 1 d 2 u . d 2 u . l2/ „ 2 


• + 


dr 2 r 2 d(p 2 dz 


+ ^ + k ^ + ju)u = 0. 


( 2 . 8 ) 


Next, we define the Q operator as 
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„ . 2 i id 2 a 2 ] 

Q- n +M+T2 ”TTTr + TT • 
A k Q [r d<t> dz 


In general, the azimuthal coupling is ignored in a Nx2D model by neglecting the —r 

a^ 

term. The influence of azimuthal coupling may be expected to be of second order, so the 
neglect of such influences here is consistent with the treatment of small perturbations. 


The wave equation now becomes 


i+k 2 Q 2 u=0. 


( 2 . 10 ) 


Defining u in terms of its incoming and outgoing components, Tappert (1977) showed 
that we may write it in the form 


=-^=(u + +u_) , and ^~ = ik 0 JQ(u + -u_) . 

V<2 dr 


( 2 . 11 ) 


The wave equations defining the evolutions of the outgoing and incoming fields can then 
be shown to satisfy 


a« + _ 1 dQ 

-*-~ik 0 Qu + +-—u_ , and 
dr 2 or 


a u. 1 dQ 


( 2 . 12 ) 


This shows how the fields form a set of coupled, one-way parabolic wave equations. 

Since the employment of an operator, Q, is not algebraic and therefore 
inconvenient, we invoke the wide-angle parabolic approximation (Thomson and 
Chapman, 1983) and define Qwape, 


Q~Qwape ~ J 1+ j2^2 + ( n2+ / / H _1 ’ 


(2.13) 


We can re-write it in terms of a Hamiltonian differential operator, H v 


Qwape ~ 1 (^ + — 1 H WAPE (r,Z, - 2 ) » 


(2.14a) 


8 



where, 


H 


WAPE 


12 \ 


r,z. 


dz 2 


= T + U 


is written in terms of kinetic and potential operators. 


T = l- 1+ 


1 d 2 
kl dz 2 


and 


U = - 


(n 2 + juY 2 ~ 1 


Since we expect n to be small, 


2 n 2 n 


(2.14b) 


(2.14c) 


(2.14d) 


(2.15) 


where we have also defined 


U(r,z) = -(n(r,z)~ 1) 

as the unperturbed potential function. Equation (2.14a) then becomes 

d 2 . 1 ju(r,z) 


Qwape ~ 1 HwAPE(r,Z, dz 2)+ 2n(r,z ) 


(2.16) 


(2.17) 


since H WAPE (r,z, ^ 2 ) — . 


i? 

3z 2 


Now, —^ » ^Qwape _ ^ dft s j nce ^ strongest range dependence is assumed in 
dr dr 2h dr 


/i. Then, the evolution of the outgoing and incoming fields from Eq. (2.12) becomes 

(2.18a) 

(2.18b) 


du. .. ., f. , 1 dju , 

i^ 0 m + — ik 0 H WAPE u + + ^ ^ u_ , and 


. a 1 

— ~ -ik 0 u_ + ik 0 H WAPE u_ + m + . 
dr 4n dr 

From Eq. (2.18a) and (2.18b), the two fields are seen to couple through the third term on 

the right-hand-side of each equation. We now invoke the multiple-forward-single- 
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backscatter approximation. This assumes that the forward propagation can only scatter 
into the backward (reverberant) propagation direction once. Since energy initially only 
propagates outward, that is, only u + exists, it is reasonable that this will contribute to u_. 

du 

This is seen in Eq. (2.18b) through the scattering term -r— . With the single- 

dr 

backscatter approximation, we then say that the secondary scatter influence on u + by 
du 

—u_ in Eq. (2.18a) is negligible. Thus, the forward propagation is assumed to be 
or 

unaffected by the small scale perturbations. The influence of is only to generate source 
terms for the backscattered field. 


1. Forward Propagation 


We can now obtain u + by defining u + = y/e ik ° r so that Eq. (2.18a) (ignoring the 
du 

—u_ term) becomes 
dr 


duf 

+ ik 0 H WAPE y/ = 0, 


(2.19) 


which defines the evolution of the forward propagating PE field function if/. The solution 
to this equation is obtained from the MMPE model with the long wavelength fluctuations 
built into n 2 . 


2. Backward Propagation 

In similar fashion to the forward propagation problem, we define an incoming 
field function, u_ = <t>e~‘ k ° r , so that Eq. (2.18b) now becomes 

= -ik 0 < D+, ( 2 - 20a > 

dr 4h dr 

noting that the second term on the right-hand side is the forward coupling term previously 
derived. The forward coupling term /, is identified as the inhomogeneous term in our 
solution for <I>, 
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(2.20b) 


J - ^ jfyf vfr ik< ‘ u ' Jrr) 

4 nd/ ¥ 

Recall that y/ and d> are the field functions of the two-way propagation with the 
long wavelength components built-in. The inhomogeneous term, I, then defines the 
coupling from scattering due to the short wavelength perturbations. In the next 
subsection, we will introduce the point source and then cast the incoming and outgoing 
fields into the form of the sonar equation and see the two-way coupling contributions 
likened to transmission losses but with a wave scattering strength. 

3. Forcing Function 

Up till now, we have only dealt with the homogeneous wave equation 

V 2 p + k 2 p = 0, (2.21a) 

with no hint of a source function. If a point source is now introduced, Eq. (2.21a) 
becomes 

V 2 p + k 2 p = -47tR^P Q 8(r -r s ), (2.21b) 

where P 0 is the source pressure at a reference distance . This may be expressed in 
terms of the Green’s function, g(r,r s ) , which is the transfer function from a point source 
to a point receiver defined by (Jensen, et al., 1994) 

(y 2 +k 2 )g(r,r s ) = -S(r-r s ), (2.22a) 

such that the pressure is 

p(r, r s ) = 47tR 0 P 0 g(r, r s ). (2.22b) 

Extending this argument to an extended source, we can apply the idea of superposition of 
a number of point sources to say that 

(V 2 + k 2 )p(r ) = f(r) = -47t^w(r s )8(r - r s )d 3 r s , (2.23a) 

then 
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p(r s ) = 4^1 w(r s )g(r,r s )J 3 r s . 


(2.23b) 


We can also define the Green’s function alternately in terms of the PE field 
function yrby letting 


4^g(r,F s ) = J%(r,r s )e ik r 
\ r 


(2.24) 


The corresponding inhomogeneous PE field function (Eq. (2.19)) for the outward 
propagating field due to a point source becomes 


^ + ik 0 H WAPE ¥ = ~4x^e*° r S(r - r s ) 


(2.25) 


since \j/\ = 1 is normalized at r = R 0 . Similarly, the inhomogeneous form for the 


scattered (or backward propagating) field is 


3 $ 

— -ik 0 H WAPE ® = -473. 


(2.26a) 


where the inhomogeneous source term is due to the scattering of the forward propagating 
field, defined by 


3 =- —^~y/e ,ko(r+r ' ) , 

16 Tin dr 

which acts like a point source function locally at the scattering point. 


(2.26b) 


We now need to define some geometrical setup to aid in the understanding and 
development of the effects of scatter. Defining the location of the scattering point as r s 

and the volume of the scattering region as V s , then the received scattered (backward 
propagating) field at location r R will be 


P-(r R ) = 4^j3(^)J ._^_. e'^~ r5l g(4,r y >fV 5 . 


(2.27) 
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Since the Green’s function is the transfer function between two given points, we can say 
that g(r R ,r s ) is the pressure field Green’s function from each point within the scattering 

region to the point receiver location at r s . By reciprocity, we say that 

g(r R ,r s )-g(r s ,r R ). ( 2 - 28a ) 

Also 


then 


g(r s ,r R ) 


1 

4/zKo 




(2.28b) 


P-(r R ) = - 


1 

•1 dyU 

*0 

16^ J 

s” d/ \ 

r s ~r R 


(r r ,F D W^'% 3 n 


(2.29) 


1 


,, -j-T^P + rP + R d ' r s > 
16^? 0 * n dr 


where p +T is the forward propagated field from the transmitter to the scattering point and 
p +R is the forward propagated field from the receiver to the scattering point which by 
reciprocity is the same as p _ R . 


1 dll 

We now examine the integrand term ——7 . Since h 2 is defined in terms of the 

n or 


step function, 


n 2 =nl(l-H(z-z b -n l ))+n 2 b H(z-z b -T)i) , 


(2.30) 


then with H(0) = —, the value of n 2 at the interface is just the average value, 
2 


n 2 (z = z b +r} l )=\(nl+n 2 b ). 


(2.31a) 


Solving for n gives 
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(nl+n 2 ) 


(2.31b) 


0 * ( n l ~ n l -Zb-Vi), 


(2.32) 


where we have ignored the slower range variations in z b , tj,, nl and n 2 b over the 
scattering region. Substituting the expressions in Eq. (2.31b) and (2.32) into Eq. (2.29), 


we obtain 




(2.33) 


_ V2 n*-n 6 2 f pTk _ n . 2 , 

Ifirc/?. M P™P+Rb a r S . 


where p +Tb and /? +Ri) are the values evaluated at the bottom interface. Note that the range 
derivative of tj s is in the direction from the scattering patch to the receiver. 

4. Statistical Treatment of the Scattered Field 

Of interest is the mean squared reverberation from different roughness 
realizations, which from Eq. (2.33) is defined as 


2_ \n w ~n b 


-JJ( •> vm 


(2.34) 


V'-' / (16^) 2 nl + n 2 b JJ\|3r 

where r s and r s indicate two different boundary realizations. Consider the statistical 
characteristics of the roughness. We note that the small scale roughness must have zero 
mean, i.e. (tj s ) = 0 and has an amplitude spectrum of S %ZD (K,L) for a given realization, 
such that 


V s (*, y) = J J S*, 2 d (K, L)e iKx e iLy dKdL , 


(2.35) 
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where K and L are the horizontal wavenumbers in the x- and ^-directions, respectively. 
Here, the amplitude spectrum, S ns lD (K,L), refers to the 2-D roughness of only the 

small-scale perturbations, therefore the subscripts ‘r/ s ,2D\ The rms roughness is then 


given by 


rj s | 2 } = JJ dKdL\\ dK'dL' S^ 2D (K, L)S^ 2D (K\ 


Assuming a wide-sense-stationary spectrum in space, then 


(2.36) 


(|>7.f} - JJ dKdLW naD (K,L), (2 ' 37) 

where W Vs2D (K,L)-^S %2D (K,L)^ is the power spectrum. The spatial correlation 


function for the roughness can then be defined as 


<17,(*, yH<.x-x',y- ft) = ||dKdL W, i 2d (K, L)e xi ~ J) e u ’-» 


(2.38) 


Hence for a given realization, 


= II klW VtM e ,K ^%^dKdL, 


where 


(2.39a) 


k r =K cos^ , (2.39b) 

and </) is angle in the direction of r towards the receiver relative to the (x,y) coordinate 
system. The following diagram shows the geometry being developed here. 



Figure 2.1- Scatter Geometry 






Conversely, the alternate-domain realization may be expressed as. 


IK 


dV s (x,y) dt] s (x-x,y-y) 


dr 


dr 


)e~ iK{x - x 'V iL(y - y ' ] d{x - x)d(y - y) 


(2.39c) 


= *X.2 d(K’D 


The above equations may be considered as the transform pairs for achieving a bottom 
realization. 


Let us now set up a spatial representation of the notations that will be used in the 
ensuing arguments. Consider two scatter points, r s and r s located at their respective 

distances from the transmitter, T, and the receiver, R. These are shown in the following 
diagram. 



The backward propagation or scattered field in Eq. (2.34) may now be written as 




drj s (x,y) dt] s (x,y) 

1 1 

dr R dr R 

4 r T r R A l r T r R 


—— Wjb ( r T )y / Rb ( r R ( r T Rb ( r R ) (2.40a) 




where 



(2.40b) 


16 




We have thus far assumed that the fields vary slowly in range, and over a given 
scattering patch the *— amplitude terms also vary slowly. Furthermore, since rj s and 

y r T r R 

decorrelate rapidly for r T £ r T and r R £ r ' R , then we may assume that these are the 
dr 

same in the amplitude terms and the locations of evaluation. Hence, the field functions 
become independent of the integral over the scattering patch so that 




(2.41) 


jy< 


drj s (x,y) dT] s (x,y') 


dr R dr R 


e *°(*+*-*-<> ]d 2 r s d 2 r' 


where r T = —{r T + 4), r R = —(r R + r R ) are the mean values of the distances from the 
2 2 

scattering patch to the receiver and transmitter. 


Introducing difference variables, r" = r r ' - r T and r R =r R -r R , and assuming the 
roughness is isotropic, W Vs 2D (K,L ) = W Vs 2D (k ), we can write the above integral as 




drj s (x,y) drj s (x-x',y-y') 


dr R dr R 


-^r^ d 2 r gd 2 r ; = J k 2 W %2D (k)d 2 r s (2.42) 


where k = & 0 2 [(cos (j) T + cos <t> R ) 2 + (sin^z> r + sin) 2 J. We now assume that over the small 
scattering patch, k 2 and W Vs 2D (k) vary slowly so that 


< 2 ' 43) 

T R 

and A \ is the area of the scattering patch on the bottom interface. The new geometry 
based on the mean distances from the scattering patch is shown in the following figure. 
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Figure 2.3 - Scatter Geometry Based on Mean Distances from Scatter Patch 


To simplify the analysis, we realign the diagram to have the receiver line coincide 

with the x-axis, i.e. (f> R = 0. Then k r = k = k 0 cos</> T and k = 2£ 0 cos-^-. Note that this 

wavenumber indicates that the spectrum is evaluated at the Bragg wavenumber. The 
following figure shows the re-oriented scatter geometry. 


i 

Receiver, r R 

V 

Transmitter, r ? 

lid Scatter patch 

A \ 


Figure 2.4 - Re-Oriented Geometry for <p R = 0 


For monostatic scatter then, we have 




(2.44) 
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where r = r T =r R and \y/ b \ = \y/ n | = \y/ Rb 


Now writing 


f 

n b =n w +An = n w 

\ 



\ 

) 


- n w (\ + 8) where 8 « 1, then 



(2.45a) 




= n 


* 2nlS 2 (l -S)~ 2n 2 S 2 = 2(A n ) 2 . 


2 ( 1 +< 5 ) 

_ , . c n c 0 CoCc^-cJ _ c 0 Ac 

Furthermore, A n = n b -n w = -=-= , 


so 


c b c„ 


C Cu 

wb 


C Cu 
wb 


tuf- 


C Cu 
wb 


'Ac } 2 

\ c ° J 


, since c w ~c fe ~c 0 . 


Finally for the monostatic case, the backscatter field becomes 


v c °; 


(2.45b) 


(2.46) 


5. Mean Squared Reverberation Pressure Level, RPL 

Of interest is the mean squared reverberation pressure level, RPL. Reciprocity of 
the one-way fields has been invoked in order to write the solution in terms of the 
outgoing field. Using the solution for the incoming field in Eq. (2.46), we can now cast 
this into the form of a SONAR equation. We do this by considering the contributions of 
the source level, directional characteristics of the source and receiver and the area of 
ensonification being captured in that region of the beams. Hence assuming that the field 
functions do not vary rapidly over the scattering region and the spectrum of the roughness 
is isotropic, the monostatic reverberation pressure level can be expressed as 


( 

RPL = 10 log 

V 



= SL + DI t + DI R +10 log 




RLfr'V ’ 


(2.47) 
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where SL is the source level, DI T/R are the directivity index for the transmitter/receiver 
respectively, and A\ is the ensonified area of the bottom interface scattering. The 

reverberation loss has also been correspondingly defined as a general function of either 
the boundary or the volume. This reverberation loss analysis applied to the boundary, 
from the analysis in the previous section, is defined by 


RL b =- lOlog 


M 2 




V u J 


(krH.ww 


= -101og 


7(kr 


-10 log 


2 


(8 ^ 


W ns , 2D (2k 0 ) 


V u J 


If the field has fully saturated statistics, which is often true, then 

(w 4 Hw 2 ) 2 . 


Then, 


RL b = -20 log 


>4 


-10 log 


' 2 ^A 




W Vs , 2D (2k 0 ) 


v u j 


= 274-10 log 


kl (Ac x 


yin 


v c °; 


W ns , w (2k 0 ) 


(2.48a) 


(2.48b) 


(2.49) 


where TL b is the transmission loss from monostatic source/receiver to the scattering 
patch at the bottom. Finally, we can write. 


RL b = 2TL b - S b , S b =10 log 


^ Ac ''' 2 


32 n 2 


\ c ° J 


K.™( 2k o) 


(2.50) 


where S b is the full-wave scattering strength due to the small-scale interface roughness, 
and W v ^ 2D (lk 0 ) is the 2-D spectrum of the interface roughness evaluated at the Bragg 
wavenumber. 
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One should note the following points: 


• If the statistics of the interface are homogeneous over a large portion of the 
area of interest, the interface reverberation structure is controlled by two-way 
transmission loss. 

• Because the two-way transmission loss is dominated by forward and 
backward propagation over long-wavelength features, the statistics of this 
signal should be related to the low spatial wavenumber-space statistics. 

The statistical model to be used for the interface roughness will be treated in the 
next chapter. 

B. VOLUME SCATTERING 

Having trudged through the analysis of the bottom interface extensively, the 
volume scattering development will take on the same general form with the exception 
that the index of refraction varies three-dimensionally in space. We will specifically 
consider for now only the theoretical development of sound speed variations and the 
overall effects on scattering and reverberation. The other element affecting volume 
scattering is density variations in the volume, which will not be treated in this thesis but 
will be left for follow-on work. 

We start with the wave equation in Eq. (2.1) and the index of refraction 
perturbations in the spatial scale in Eq. (2.3), 

V 2 p + £ 0 2 n 2 p = 0, (2.51) 

and 

n 2 =n 2 w a-H(z-z b )) + n 2 b H(z-z b ), ( 2 - 52 ) 

as before, but now the bottom (volume) index of refraction varies three-dimensionally in 
space, n b (x,y,z ). As before in the statistical treatment of the interface spatial variation, 
sound speed variations may be defined by (Yamamoto, 1995) 
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C b (x, y, z ) = c b0 (1 + bz + 8(x, y, z )), 


( 2 . 53 ) 


where c b0 is the mean bottom sound speed at the interface, b = — is the normalized 

C b0 

gradient of the bottom sound speed, and 8 is the zero-mean random perturbation, i.e. 
(<?} = 0 and 

(8(x,y,z),8(x',y',z)) = JJJW^ 30 (K,L,M)e iK{x ~ x) e iUy ~ y) e iM< - z ~ z) dKdLdM , (2-54) 

where W S 3D (K,L,M) represents the 3-D power spectrum of the fluctuations. 


Again we break this up into long and short wavelength components by having the 
random perturbations split as 


s=s,+s,, 

such that Eq. (2.53) now becomes 

c b (x, y , z) = c b0 (1 + bz + S ,) + c b0 8 s 
= c b (x,y,z) + c b0 S s (x,y,z). 

This can be incorporated into the index of refraction by writing 




n b = 


\. Cb J 


c b0 (1 + bz + S,+ 8 S ) 2 


-^a-2 S s ) = n 2 b+ M b , 


where now 


(2.55) 


(2.56) 


(2.57a) 


ju b =~2n;S s (2.57b) 

and hi is, as before, the description of the environment including all long wavelenth 
component perturbations. 

The general description of the environment of Eq. (2.52) is now 


h 2 = <(!■- H(z- z b -r)i )) + n 2 b H(z- z b -r t } t ), 


(2.58) 
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and the effects of interface scatter and volume scatter may be considered separately. Note 
that the long-wavelength contribution to the interface roughness has been included as r\ l 
into the perturbation, which in turn affects the effective index of refraction. The 
derivation now follows exactly as before arriving at the final result for the bottom 
interface in Eq. (2.33), 


P-=- 


f-i 


1 dju b 


P+tP+r d r s > 


(2.59) 


16^ n dr 

where the volume integration occurs over the sediment volume scattering region, i.e. 

Z>Z b +TJ l . 


We now have 


^=*-2 n 2 


dr 


dr 


which remains a volumetric quantity, so 


(2.60) 


f - j 3 

P-=T7-^) n -r:P+TP+Rd h 


\67tR Q l dr 


(2.61) 


Taking the same approach as the interface in Eq. (2.34), we see that the mean squared 
volume reverberation is then 






(2.62) 


P+TbP+RbP+TbP+Rb)d r $ d Tj . 

As before, we write the pressure fields in terms of the field function, if/ (see Eq. (2.40a)), 
and argue that they do not change significantly in range within the scattering region. 
They do, however, vary considerably in depth, so 


(\p.f) = . 1 v j-j-Jj (w T ( r r . z)¥ T(r T , z)Wr (r R , z)w\ (r K , z)n(r s )n(r' )} (2.63a) 

' ' (o^K 0 j r T r R zz 


X 


11 


xyxy 


, \ dSc . r t t \| 

dr dr 


e ilc 0 (r T+rB -4-rl) 
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Recall that r T =^-(r T + r') and r R = —(r R +r' R ). The inner term is evaluated over 

A 

difference variables (as before in Eq. (2.42)) r” and r" to give 


II 




■ V ’ ^ir~ (X ~ X " ,y ~ y "’ z ~ Z ^ 


)e' iko(r ^ ) )d 2 r s d 2 n 


(2.63b) 


= jk*W Ss ' 2DH (K',L")d 2 r s 
where K", L” and k r are as before, but 

W, sj2DB (£,L)^ 3D (tf,^ (2.63 c) 

is the 2D spectrum of the small scale perturbations in the horizontal dimension. 


If we now assume that the spectral density is essentially a delta-correlated 
function in depth and horizontally isotropic, then we can say that 
w s s ,2D„ (K,L) = W Ss2Dh ( k) = w s s ,3 D (k,0)S(z-z) , and neglect variations in W Ss2Dfi over 
the scattering region. Then, assuming a monostatic geometry, 

(I /? ‘I 2 ) = (8^j 7‘ f (” 2 ( r ’4^4W r ’ z )| 2 ) (2 * 64) 

x k?W Ss 2Dfi (k)d Z AA , 

AA being the horizontal area extent of the scattering volume, much akin to the scattering 
patch for the interface reverberation development. 

For backscatter then, 

(H 2 }=(—^ 7 - / [n\r,ziwAr,zf\y/ R {r,z)\ 2 ) (2.65a) 

xk 2 0 W Ss2Dii (2k 0 )dzAA. 

If we now assume that the spectrum W Ss2Dh does not change significantly with depth, 
then 
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(2.65b) 


(k| 2 )4^;l W Ss ' 2Dh (2Jc 0 )AA ](n 2 (r,z)\Hr,z)\yz. 

X \ J Zb+T h 

Finally, we will assume that ^n 2 (r,z)|^(r,z)| 4 ^«2^n(r,z)|^(r,z)| 2 ^ , similar to our 
previous treatment of fully saturated statistics. 


For the volume reverberation problem, we see that we cannot write the RPL 
simply in terms of transmission loss (TL) but instead must integrate over depth the 

quantity ^n(r,z)|^(r,z)| 2 ^ at each range, r. We shall then use a form of ergodicity and 

write 


J^n(r,z)|^(r,z)| 2 ^ dz » n(r, z)|^(r, zfdzj , (2.65c) 

and since n and y/ are based only on the long wavelength components, we will simply 
say that we can use {^h{r,z'jfl/{r,z'j[ dz) . Then as in Eq. (2.47), we may write 


PPL = 10 log 




_2 
Pref 

V ' ' J 


= SL + DI t +DI r +10 log 






-RL V , 


(2.66a) 


where 


RL V =- 20 log 


■p °°\h(r,z)\v(r,z)fdz 


z>z b +v. 


-10 log 


k 2 

ZQ w 

2 S s ,2D„ 


yin 


w 


(2.66b) 


In the above expression, W S} 2Dh (2k 0 ) is the 2-D horizontal spectrum of the 

volume fluctuations, which is assumed isotropic and independent of depth (as in the 
interface reverberation problem) evaluated at the Bragg wavenumber, and the depth 
integral extends from the bottom interface downwards to negative infinity. Note that 
volume reverberation loss cannot be expressed simply in terms of the two-way 
transmission loss as before. 
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One should also note the following few points: 

• The reverberation loss is not independently influenced by the forward 
propagating fields, but is a magnitude scaling or weighting function of the 
local index of refraction (positive definite). 

• If the statistics of the volume are homogeneous over a large portion of the area 
of interest, the volume reverberation structure is controlled by this non-linear, 
two-way transmission loss between sound speed structure and acoustic field. 

• Because this weighted two-way transmission loss is dominated by forward 
and reciprocal propagation over long-wavelength features, the statistics of this 
signal should be related to the low spatial wavenumber-space statistics, as in 
the case of interface reverberation. We will however expect to see smoother 
changes (not abrupt) as compared to the interface reverberation statistics. This 
is because the statistics apply non-linearly here providing only the mean- 
square positive-definite weighted average of the acoustic interactions with the 
environment. 

The statistical model to be used for the volume sound speed variations will be 
treated in Chapter HI. 

C. TIME-DOMAIN ANALYSIS 

1. Time-Domain Analysis of the Interface 

Up until now, the focus has been on CW analysis. The previous analyses for both 
the interface and volume have suggested that the structure of the reverberation is due 
mainly to the effects of propagation. To predict the effects of a pulsed signal instead of 
the CW signal, the MMPE/MMPEREV model is run over a spectrum of frequencies. 
Subsequently, Fourier analysis will give the time-domain structure of the pulse 
propagation in time. 
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In shallow water, propagation of a signal from a point source will experience 
multipath effects. In the time-domain, a scatter patch receiving multipath arrivals will 
scatter the received signal at different arrival times. The scattering patch, by reciprocity, 
scatters the entire received signal back through the same multipaths (monostatic 
geometry assumed) that each signal has traversed, leading to the receiver receiving 
multiple-multipath signals at different times. Thus, the scattering patch is a source, 
producing time-delayed replicas (with proper amplitude weightings) back to the receiver 
through the multipath structure of the water column. The geometry of the development 
may be represented by the following diagram. 



As the source transmits a signal out to the scattering patch, the scattering patch 
receives signals along all the multipaths (in the figure above, just three multipaths). The 
scattering patch scatters each signal and, by reciprocity, these scattered signals follow the 
same (three) multipaths back to the source location. Recall that the Green’s function is 
the transfer function between the source and the receiver and it does not matter which 
direction the propagating wave is travelling. Hence, from a signal point of view, these 
may be represented as depicted in Figure 2.6. The total reverberation signal from this 
scattering patch would then be the coherent sum of the three lower panels. 
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Jigrtal from Path 1 


Path 2 



In general, if the receiver and source are not co-located, it may be useful to 
consider the arrival times of the same signal from the source and receiver to the same 
scatter point/patch. By reciprocity, we may say that the receiver has “propagated” the 
signal towards the scatter patch. The geometry has been adapted from Smith and 
Cushman (1997) as shown, in Figure 2.7. 



Figure 2.7 - Two-way return from a Scattering Patch 
These two propagations may be considered independent and in order to compute the total 
reverberant field at the receiver, we may convolve the two field functions as explained by 
Smith and Cushman (1997). With an understanding of the pulsed propagation and 













multipaths, the one-way time-domain pressure due to a pulsed signal at the scattering 
patch is 

P+T/Rb( r m’ty = J 

where the forward propagated pressure field from the transmitter/receiver at the interface, 
p +T/Rb , is specified at range r m , which is the range at the m h range step. The field- 

function evaluated at the interface at range r m and at frequency / may be extracted from 
the forward propagation of the broadband field computed by MMPE. Notice that in order 
to arrive at the required pressure in the time domain, we have to integrate over all the 
frequencies, and the above expression is really the Fourier transform of the field in the 
frequency domain. 

Extending the argument further, the travel time of the two-way pressure field at 
the receiver is now the convolution (Smith and Cushman, 1997) of two, one-way fields in 
the time-domain propagated forward. This is 

P 2 -way,,b C T m ’ 7 ) = J P+Tb ( r m > *)P+Rb ( r m ’ 7 _ 0^ * (2.68a) 

where we recall from Eq. (2.33) that p +Tb and p +Rb are the forward propagated fields 

from the transmitter and receiver to the scattering point evaluated at the bottom interface, 
respectively. Note that the forward propagated field from the receiver to the scattering 
point is by reciprocity the same as p_ Rb . We also note that the transmitter and receiver 

need not be co-located in depth. 

Applying the above approach, we see that the time-convolution of the two field 
functions is also the frequency multiplication of these functions in the frequency domain. 
Hence, when the forward propagating fields are already available in the frequency 
domain, we may then compute the two-way field in the frequency domain from the 
interface as 

P 2 -way,b ( r m . /) = P+Tb <T m > f)P+Rb ( r m ’ f) ’ (2.68b) 


4Z 


y^T/Rb ( r m’f) e 


.27? 

:——r 
Co 


I e- iW df , 


(2.67) 
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where 


P + n ( r m ./) = ~j=¥ + n ( r m ,/> iV , and 
Jr„ 


(2.68c) 


P+Rb( r m’f^~ I - V^+Rb( r m’f) e ° » 

v r - 

at a particular bottom depth, z b {r m ) , at range r m , and at each frequency. By applying Eq. 


(2.67), we can then arrive at the two-way travel time structure of the reverberation loss 
(RLb) due to a single bottom patch, i.e. 


P.(r m , 1) = Aj p 2 ., vt (.r m ,f)e-‘ 2 *df, < 2 ' 68d > 

where the constant A has been included to incorporate all the other factors needed to 
define reverberation loss. 


We then continue our propagation through the entire water space of interest and 
arrive at the reverberant field at each range step. To see the total field at the receiver, we 
coherently sum up all the pressure values from different range segments, r m , by matching 
up the discrete arrival times, t n , 

P-W-l>.<r..O <269) 

m -1 

where p_ is the total pressure at the receiver received at time t„. Note that this matching 
requires a coordinated scaling between range step size and time step size, as will be 
explained in Chapter IV. 


2. Time-Domain Analysis of the Volume 

The volume analysis is a 2-D problem in that the reverberation is not directly 
proportional to the two-way transmission loss, but requires proper weights through the 
entire volume when computing the reverberant field at the receiver. The development 
will follow closely that of the previous sub-section on the interface. 

Firstly, the reverberant field due to each depth/range point is computed using the 

frequency-multiplication of the source-patch and patch-receiver (by reciprocity, the 
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forward propagation used) propagating fields, and then applying a Fourier transform to 
the time-domain. This is 

Pi- way (r m ,Z,f) = n{r m , z)p +T (r m , z, f)p +R (r m , z, /), (2-70a) 

where now the two-way reverberation signal is computed for every grid point of interest 
(z> z b always) at a particular frequency, /, and n(r m ,z) is the local index of refraction 
at the grid point. The justification for multiplying by the local index of refraction is to 
provide the same weighting used in the CW analysis. Fourier transforming gives 

Pi- way ( r m = jp 2 _ way (r m , z, f)e~ a ^df . (2J0b) 

The 2-D problem may now be reduced to a 1-D problem by performing a vertical 
integration to arrive at a single reverberation time-series, i.e. 

CO 

P-(r m ,t) = B jp 2 _ way (r m ,t)dz, (2.71) 

J>2 t 

which is the coherent sum over all depths below the interface at range step m. Note that 
the constant B has been introduced to account for all remaining terms in the reverberation 
loss, RLv. We have now ‘collapsed’ the volume in 2-D to 1-D, giving a single set of time- 
series which can then be matched and summed, 

(2 - 72) 

m=\ 

giving a single time-series in p_(t„) which is the two-way time-domain field defining the 
volume reverberation loss at the receiver at time t„ due to the entire volume in the space 
of interest. 

In the post-processing analysis to be discussed in Chapter IV, we will implement 
the above discussions for both the interface and volume to arrive at the time-domain 
analysis. 
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III. GENERATING REALIZATIONS - IMPLEMENTATION OF 
THE REVERBERATION PROBLEM IN MMPE 


Having developed the theoretical basis for modeling the interface and bottom 
variabilities, we now have to define statistical methods or models to generate the 
interface and sound speed variations. We first developed the theoretical aspects of the 
models used for generating both the interface and volume perturbations, then 
implemented the algorithms in MATLAB to enable checking of the theoretical 
developments before integrating them as functions into the MMPE model (Smith, 2000) 
written in FORTRAN 95. Much of the theoretical development for generating 
perturbations to the interface is based on the work of Goff and Jordon (1988) while the 
model for the volume is based on Yamamoto’s (1995) work. 

We should note that generating good realizations would form an important first 
step in the effort to model shallow water reverberation and their effects on sound 
propagation. Of particular importance are the “tweaking” that must be afforded in root- 
mean-square values when generating these realizations so that we are able to study the 
effects of different bottom undulations and sediment types. 

A. SPECTRAL DESCRIPTIONS OF BOTTOM VARIABILITY 

1. Interface Roughness (Spatial Perturbation) 

From Chapter II, we have seen that the interface effects can be separated into both 
long- and short-wavelength fluctuations. We have principally neglected the short 
wavelength components in the forward propagation problem. Since the interface 
roughness is essentially 2-D, we can assume a 2-D interface spectrum of the form 

-— 5/ . *,=VF+7F, (3-D 

(l+O, 2 ) /! 
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where L corr is a correlation length scale, is a spectral exponent, and k r is the 


horizontal spatial wavenumber vector. Note that for k r » 



w VS ' 2 D (k r »L corr )~k;‘ 


The normalization factor ju is defined in terms of the root-mean-square (rms) roughness 
a 2 by requiring 


27r jW r]s2D (k r )k r dk r =(J 2 , (3.3) 

o 

which leads to 

<3Aa) 

n 2 J 

Note that k r is the horizontal component of the wavenumber, and the 2-D spectrum, 
W Vs 2D , is assumed isotropic (independent of direction). 


For this description, the exponent P n is constrained by 

2<j3 v <4. (3.4b) 

Consistent with previous reverberation work by Goff & Jordan (1988) on a rough 
interface, we shall assume a fractal character with /? = 3.5. 

For the scattering amplitude, we may simply evaluate W Vs2D {k r ) at the Bragg 
wavenumber, k r = 2k 0 . However, for the long-wavelength interface roughness, we need 
the 1-D spectrum along the x-axis, W Vs W (K ). We do this by taking the 1-D transform of 
W % 2d (K, L) along the slice y = 0, 

W VS ' 1D (K)= ]w % ' 2D (K,L)dL. (3.5a) 
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In cylindrical coordinates, we have 


where 


'W*>= 2 | -vhf W 2 (K)dk r 


A. 1 


= y<T% orr (\ + L 2 corr K 2 ) 2 2 . 


y = 


(B n \ 

- I 

{A, i) 

2 11 

1 11 

1 2 111 




m 


v 2 ; 


(3.5b) 


(3.5c) 


We now have to transform a random realization of the corresponding 1-D 
amplitude spectrum which has been scaled by a random amplitude and phase. In other 
words, we define the realization as 


*K*)=J^ iU> (jry s, <flr (36a) 

—oo 

where 

S l (K) = [w VS ' W (K)Y 2 A(K)e i6(K) , (3.6b) 

and A and 0 are random numbers for all values of K. Since the complex amplitude of 
each wavenumber component, Ae ' 9 , should exhibit a normal distribution, the random 
phase of each component can be obtained from 

6 - 2m x , (3.7a) 

where r\ is a uniformly distributed random variable in the interval [0,1]. The magnitude 
A, however, exhibits a Rayleigh distribution. Because we are in fact generating a 
realization of the power spectrum, we must consider the magnitude-squared, which has a 
negative exponential distribution. The random amplitude of each wavenumber 

component is then obtained by 
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A = V- ln(r 2 ), (3-7b) 

where r 2 is another independent, uniformly distributed random variable in the interval 

[ 0 , 1 ]. 


In practice, we could simply use 


A, i 


f w »ro=(i+4«’)' ! 1 

T}(x) 

and rescale the result by its rms value, i.e. -—-——y , or we could simply say 


(3.8a) 


(*7 2 C*)) 


W. 




m~(i 


+ L; 


2 

'corr 


k 2 ) J i 


i+± 


(3.8b) 


where L corr and /3 n had to be specified. This was the generic spectral model used in 

generating the realizations first in MATLAB and then for implementation in the MMPE 
model. 


2. Volume Sound Speed Fluctuations 

Unlike spatial variations of the interface, volumetric variations are generally 3-D 
in nature. We will now treat the generation of volumetric variations through spectral 
realizations as in the interface. The sediment volume sound speed perturbations may be 
modeled by a 3-D volume spectra (Yamamoto, 1995) given by 

W s ^ D {K,L,M)=^y-[^{K 2 + L 2 )+M 2 )^ ( 3 - 9 ) 

where B is the spectral strength constant, (3 S is the spectral exponent, and A = — = — is 

a \ 

the horizontal-to-vertical aspect ratio describing the anisotropy of fluctuations in the 
sediment. We shall assume the major axis corresponding to M is aligned with the vertical, 
thus we neglect any “dip” angle in the structure. M is thus the vertical wavenumber. 
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From Yamamoto’s findings, we shall assume the following values: 


A~2, 

B ~ 5xl0 -4 , (310) 

A ~ 5. 

This should produce a relative velocity variance of ~10' 4 . Also note in this form that we 
really assume K = \k x where fy, =lm is simply a reference length so that K, L, M are 

unitless and units of W^ 3D do not depend on ($. 


Recalling Eq. (2.66b), in order to evaluate the scattering due to the volume 
perturbations, we need a representation for the volume spectra in the 2-D horizontal 
plane. Hence for our reverberation computations, we have 


OO 


(3.11) 


Note that we are effectively evaluating the 3-D spectra at a single depth, since it is delta- 
function correlated. Substituting Eq. (3.9) into (3.11), we obtain, 


OO 

W Ss ' 2 Dh (K,L)= jw Ss ' 3 D (K,L,M)dM 


(3.12) 


= Ml ffr (k 2 + L 2 )+ M 2 \ 2* _1 dM 

TT " 


n 


rfA+il 

m te*£\rfr> + e )l 2 2 

2 n 


2 2 —X- 

{‘t 


+ 1 


J 


For P s -2, Eq. (3.12) can also be shown to reduce simply to 


W s>m (K,L)=—]\A 2 (K 1 +I?)+M 2 ]~ 2 JM 

n 




(3.13) 
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For the values of B and T in Eq. (3.10) chosen, this reduces to 

W s ,.v,W.L) = a(K* + I?f I ( 3 ' 14a > 

where 

a = — = 5xl0" 5 . (3.14b) 

2A 

Note that the above provides the 2-D horizontal spectrum which is needed to 
evaluate the reverberation. However, to compute the forward propagation, we need the 2- 
D vertical spectrum in the (r,z) plane. In other words, we will now perform the 
integration cross-range at y = 0, giving 



(3.15) 


(3.16a) 


(3.16b) 


To generate 2-D vertical volume fluctuation realizations, we do essentially the 
same thing as in the case of the interface. Since we are looking at sound speed 
perturbations, we define a realization as 


c 0 S(x, z ) = JJ S Ss , 2Dv {K,M)e iKx e m dKdM 


(3.17a) 
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where 


S s „^K,M)A^s, m (K,M)Y 1 A(K,M)e mKM) . (3.17b) 

Notice that since we have treated sound speed perturbations in the volume in the vertical , 
we really are generating a series of vertical realizations at each range step. 

In similar fashion to the interface, we define the 2-D random phase variations as 
d{K,M) = 27tr l {K,M) , (3.18a) 

where r x (K,M) is now a matrix of uniformly distributed random numbers in [0,1], and 
the amplitude variations are 

A{K,M) = 4~Wr 2 {K,M)) , (3 - 18b) 

where r 2 (K,M) is another matrix of uniformly distributed random numbers in [0,1]. In 
practice, we use 

W Sim (,K,M)~{/?K 2 , (3 ' 19) 

where A and fi s are to be specified in generating the realization. The final result for the 
volume perturbation is then rescaled by the desired rms value. 


3. Spectral Filter 

We must note that both the interface and volume realizations will need to be 
filtered at low and high wavenumbers in order to reduce edge effects, which will 
introduce additional unwanted spectral components. The filter is also used to extract the 
pertinent portions of both the interface and volume spectra for creation of the long- 
wavelength realization affecting the forward propagation. 


At the low end, we can assume scales greater than 500m are known. As such. 


min 


2 n 
500 


= 0.1 


(3.20a) 
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At the high end, we should compute up to about the acoustic wavelength, so 



which will depend on our frequency of calculation. In a broadband computation, c 0 will 
be our center frequency. 

B. IMPLEMENTATION OF THE INTERFACE AND VOLUME 
SPECTRAL DESCRIPTIONS IN MMPE 

As the MMPE is well documented (Smith, 2000), a quick run-through will be 
afforded in order to understand the integration of the spectral perturbation model into the 
MMPE. However, before any implementation was possible, MATLAB was used to test 
the model and its algorithm. 

1. Description of MMPE 

The Nx2D MMPE model (Smith, 2000) is based on the split-step Fourier 
algorithm developed in the UMPE model (Smith and Tappert, 1994) written in 
FORTRAN ‘95. MMPE generates the environmental propagator, U op , in a subroutine 
called “ENVPROP” using the bottom bathymetric and sound speed data. These data 
points are interpolated to fit the depth points being computed at each range step. 
ENVPROP is subsequently called at each range to generate the propagator variations due 
to a range-dependent environment specified in the input files. In the case where there are 
perturbations, the program is ‘forced’ to think that it is in a range-dependent 
environment. 

Several input files are used to specify pertinent input parameters, most important 
of all being the file named “pefiles.inp.” In there, several input filenames are specified, 
one each for the source parameters, sound speed profile(s), bathymetric data set(s), 
bottom properties (like shear and compressional sound speeds, etc.), a sub-layer 
bathymetry, and the sub-layer properties. A total of seven input files are used for running 
MMPE simulations. 
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2. FORTRAN & MATLAB Implementation of Perturbation 
Models 

Before actual implementation of Eqs. (3.8b) and (3.19) into the MMPE’s 
FORTRAN code, MATLAB was used to code up and test the algorithm and effects of 
generating different roughness and perturbation rms values. These were proven 
successful, with the expected spatial and sound speed perturbations modeled effectively. 
A “tweaking” factor for each type of perturbation was used, namely the rms roughness 
and the rms sound speed variations, to allow observation of these perturbations. 

Subsequently, these algorithms were translated and implemented into the MMPE, 
now renamed “MMPEREV.” In MMPEREV, two additional subroutines are called 
within the subroutine “ENVPROP1” to introduce perturbations. “ZBGEN” introduces the 
spatial perturbation to the bottom bathymetry at each range step while “DCGEN” 
introduces the sound speed perturbation to the volume at each depth grid point and range 
step. ENVPROP1 is then made to think that a range-dependent case is present simply 
because the perturbations cause irregularities in depth and sound speed, causing the 
program to call ENVPROP at each range step. 

A random number generator was used in FORTRAN to provide the uniformly 
distributed random numbers in the range [0,1]. Each time the random number is required, 
a series is generated. However, each time the compiled MMPEREV program is started, 
the same seed is used by FORTRAN so that the same set of random numbers are 
invoked. Thus, the same sequence of random numbers is used for each run, allowing for 
controlled analysis of the results. This is described later. Later implementations to 
MMPEREV will allow for the generation of a random seed in order to effectively model 
truly random fluctuations to the intended data set. Some results of the controlled analysis 
are seen in the next sub-section. 

3. Results of Varying RMS Values to Perturbations 

We now present some results from the random fluctuations to bottom bathymetry 
in range by varying the root-mean-square values of the interface spatial perturbations. 
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Depth, m Depth, 


Interface rms roughness: 1m 



0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Range, km 

interface rms roughness: 5m 



0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Range,km 

Figure 3.1- Plot of Interface Roughness of Different rms Values 


Interface rms roughness: 10m 










From Figure 3.1 and Figure 3.2, increased rms roughness is seen as the root- 
mean-square values are increased. The longer wavelength perturbations also become 
more distinct. A cursory survey of shallow water areas, especially in the East China Sea 
where the ASI AEX experiment will be conducted, reveals that the long wavelength 
components are relatively small. Hence, an rms value of lm was selected. 

Effects of tweaking volumetric root-mean-squared perturbations to the sound 
speed follows. The mean bottom sediment sound speed selected was 1700 m/s, which is 
normal for a sand/mud layer. 

Sound Speed Profile Data Over Depth & Range, (m/s), Volume Perturbation: 5m/s 

1710 

1700 

1690 

1680 

1670 

1660 

1650 

1640 

0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 

Range (km) 

Figure 3.3 - SSP Data with Volume Sound Speed Perturbation of 5 m/s 
and Interface Perturbation of 1 m 
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Depth (m) hrf Depth ( m ) 


Sound Speed Profile Data Over Depth & Range, (m/s), Volume Perturbation: 15m/s 




1740 


1720 


1700 


1680 


1660 


1640 


Range (km) 

igure 3.4 - SSP Data with Volume Sound Speed Perturbation of 15 m/s 
and Interface Perturbation of 1 m 


Sound Speed Profile Data Over Depth & Range, (m/s), Volume Perturbations: 25m/s 



Range (km) 


Figure 3.5 - SSP Data with Volume Sound Speed Perturbation of 25 m/s 
and Interface Perturbation of 1 m 
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Apparent from Figure 3.3 through Figure 3.5 is the increasing contrast in the sound speed 
profile for an increasing sound speed rms perturbation when comparing them 
qualitatively. While it is not likely that the sediment layer is homogeneous, as in Figure 
3.3 (rms perturbation being 5 m/s), it is not unimaginable that the sound speed profile can 
vary considerably due to inhomogeneities and mixture of different sediment types. A 
median sound speed perturbation of 15 m/s (Figure 3.4) was finally used for the follow- 
on analysis. Notice also that the interface rms roughness used was 1 m. 
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IV. POST-PROCESSING IMPLEMENTATIONS IN MATLAB AND 

SIMULATION RESULTS 


The second and final step in our modeling efforts was the development of post¬ 
processing routines to extract the relevant data from two, one-way propagation data sets 
generated by MMPEREV (see Chapters II and El) and then compute reverberation and 
perform signal analyses on the derived data. We will also see that the time-domain 
convolution of the two one-way field functions will give the time-domain reverberation 
structure received for both monostatic and bistatic setups. Routines were written to 
perform signal processing on these reverberation data sets. Specifically, this chapter will 
expand on the work done to derive the following: 

• reverberation loss for broadband signals (time-series) and CW signals (range- 
series); 

• vertical correlation and peak correlation computation; and 

• spectral analyses. 

We saw from Chapter II, Sections B and C that a scaling or normalization due to 
the wave scattering strength is needed in order to arrive at the reverberation loss. 
However, we will not attempt to introduce any scaling here, as these normalization 
factors are assumed to be constants (see Eqs (2.50) and (2.66b) for CW analyses, and 
Eqs. (2.68b), (2.68c) and (2.71) for time-domain analyses). 

Discussion of the reverberation computation cannot be made without first 
defining the geometrical and environmental setup of the modeling effort. Of interest was 
the particular configuration and environmental parameters likely to be used in ASIAEX. 
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A. MODELING ENVIRONMENT AND GEOMETRY 


Since the varying spatial properties of the broadband signal are of interest, an 
array was typically chosen in collecting field data. Each element of the vertical line array 
(VLA) was assumed to be a point source and the array was located vertically in the water 
column. The MMPE/MMPEREV model provided as an output the field function for the 
entire water column and bottom volume up to a given range. The fluctuations in the 
bottom bathymetry and volume sound speed were randomly generated by MMPEREV as 
explained in Chapter HI. 

1. Multi-static Reverberation Geometry 

In the analyses, there was a requirement to develop a geometry able to support 
both monostatic and bistatic reverberation computations/measurements. As such, a 16- 
element VLA was used in the modeling effort. Figure 4.1 shows the VLA and its 
geometry relative to a given scattering patch at horizontal range r. We see that the mean 
horizontal distance of the VLA to a small scattering patch is r, in see Eq. (2.43). 



Figure 4.1- Geometry of VLA and Scattering Patch 
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From Figure 4.1, we see that the 16-element VLA spans the water column from 20 m to 
80 m with the array elements located 4 m apart in depth. A single element at 48 m depth 
(or element number 8 from the top) was chosen as the source with all 16 elements 
receiving the reverberation. 

As there was a need to isolate source/receiver pairs, MMPEREV was set to 
generate 16 sets of output files using each array element as a reciprocal source and then 
propagating the field through the entire space with the environmental and source 
parameters specified in the next sub-section. 


2. The Environment 

The maximum propagation range was selected to be 5 km. In a shallow water 
environment, the mean bottom depth was assumed to be 100 m with no mean slope. To 
perform broadband as well as CW analysis, a center frequency of 250 Hz was chosen 
with a 250 Hz bandwidth divided into 512 discrete propagation frequencies. The 
frequency span was from 124.7554 through 374.7554 Hz computed automatically by 
MMPEREV. Seven control and environmental input files were required by 
MMPE/MMPEREV to perform the computations. In order to properly define the 
environment for running the MMPEREV model, there was a need to define the baseline 
features of frequency, depth, grid (depth and range), etc. 


In order that the time steps and range steps coincide appropriately for coherent 
addition as needed, we start by seeing that in reduced time at the m th range step, we have 


T m =t m —— , 


(4.1) 


where T m is reduced time, t m is true time at range step m, r m is the range, and c 0 is the 
reference sound speed. Hence, 


t =t+-zl. 


(4.2) 


Now, true time at r m+ 1 is 
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(4-3) 


^ ^, 4.1 rr. r + A r 
? m+1 =T+-^ = T + — - 


where A r is the range step size (not the MMPEREV computation range step size, which 
will be notated as Sr ), which leads to 


^m +1 


The frequency step is defined by 

A/= —, (4.5) 

nf -1 

where nf is the number of frequencies we would like to compute in a broadband scenario 
with BW being the bandwidth selected. The time increment then becomes 

Al = ——— = - L. (4.6) 

«/-1 (n/-l)A/ BW 

Suppose we want a bandwidth of 250 Hz, a center frequency of 250 Hz, and a 
frequency step size of about V 2 Hz, then selecting 


A t = -= 4 ms, (4.7) 

250 

we get 

nf = 512, (4.8) 

250 

A/= — = 0.49 Hz. 

511 

Note that nf has to be a radix-2 integer because of the Fast-Fourier Transform (FFT) that 
was used to compute spectral components in MMPEREV and in the MATLAB post¬ 
processing. 

We now constrain the range step to be an integer multiple of the MMPEREV 
computational range step size by defining 
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Ar = kxSr, (4.9) 

where it is an arbitrary integer to constrain the relationship with the MMPEREV 
computation range step size, Sr . By forcing this constraint, we have defined the spatial 
resolution of the generated solution to explicitly overlap the temporal resolution. This 
will allow us to easily combine solutions from different range steps coherently in the 
time-domain. 

To resolve time, we also define 


— =jxAt, 

c 0 


(4.10) 


where j is another arbitrary integer to constrain the time-range relationship for the 
purpose of computations in MMPEREV. For At = 4ms, selecting j = k= 1 and 


c 0 =1500m/s, we get 


Ar = 6m, 
Sr = 6 m. 


(4.11) 


Hence for a given maximum range of 5 km, the number of range steps to be computed is 


5000 

nr =- 

6 


« 833 


(4.12) 


in MMPEREV. Note that the maximum computational range will then be 
=4.998 km. 


Finally, we select the number of total computational depth points to be 256, which 
would also determine the FFT-size for the split-step Fourier algorithm, also a radix-2 
integer. This provides a vertical spatial resolution for computation of roughly 3 m. The 
following table summarizes the environmental parameters used in the modeling in 
MMPE, noting that much of the parameters used for the bottom and deep-bottom 
properties were assumed. 
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Filename/Parameter 

Value 

Remarks 

Main Control File: pefiles.inp 



Number of depth points 

256 

Radix-2 integer required 

Minimum depth 

0 m 


Maximum depth 



Number of range steps 

833 


Minimum range 

0 m 


Maximum range 

5.0 km 


Range step size 

6 m 


Maximum computed depth 

400 m 


Reference sound speed 

1500 m/s 





Source depths 

Varying 

Array elements at 20, 24,28, 32, 
36,40,44,48,52,56,60,64, 

68,72,76 and 80 m depths. 

Center frequency 

250 Hz 


Frequency bandwidth 

250 Hz 


No. of Frequencies 

512 

Radix-2 integer required 

Sound Speed File: pessp.inp 



Water column sound speed 

1500 m/s 

Range independent 

No. ofSSPs 

1 


Bathymetry: pebath.inp 



Mean bottom depth 

100 m 

Range independent 

No. of depth points 

1 


Bottom properties: pebotprop.inp 



Bottom sound speed 

1700 m/s 


Sound speed gradient 

0 


Relative density 

1.0 

No density variations 

Compressional attenuation 

0.2 dB/km/Hz 


Shear speed 

0 

Not modeled 

Shear attenuation 

0 

Not modeled 

Deep Bottom Bathymetry: pedbath.inp 



Depth 

3000 m 


Deep Bottom Properties: pedbotprop.inp 



Deep bottom sound speed 

2000 m/s 


Sound speed gradient 

0 


Relative density 

3.0 

No density variations 

Compressional attenuation 

0.25 dB/km/Hz 


Shear speed 

0 

Not modeled 

Shear attenuation 

0 

Not modeled 

RMS Perturbations (input to MMPEREV 
during run) 



Interface 

1 m 

As explained in Chapter III 

Volume 

15 m/s 


Table 4.1-MMPEREV] 

iiput Environmental Properties 
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3. MMPE and MMPEREV Runs 

With the propagation environment specified, MMPEREV was setup to make 16 
runs, each for a different source depth from 20-80 m at 4 m increments in depth. The 
ability to propagate the same signal out from the various receivers to the scattering patch 
is based on reciprocity as explained in Chapter II, and using Eq. (2.68b) we are able to 
arrive at the reverberant field. As such, 16 output binary files for the configuration above 
were computed and stored. In order to make comparisons with the unperturbed data, 
MMPE was also setup to make 16 runs for the same array element depths based on the 
same environmental conditions shown in Table 4.1. The following table summarizes the 
filenames used. 


Receiver 

/Source 

Depth 

MMPEREV Output File 
(perturbed data) 

MMPE Output 
File (unperturbed 
data) 

20m 

peoutzbrms 1 dcrms 15 s20.bin 

peouts20.bin 

24m 

peoutzbrms 1 dcrms 15 s24.bin 

peouts24.bin 

28m 

peoutzbrms 1 dcrms 15s28.bin 

peouts28.bin 

32m 

peoutzbrms ldcrms 15s32.bin 

peouts32.bin 

36m 

peoutzbrms Idcrms 15s36.bin 

peouts36.bin 

40m 

peoutzbrms 1 dcrms 15 s40.bin 

peouts40.bin 

44m 

peoutzbrms Idcrms 15s44.bin 

peouts44.bin 

48m 

peoutzbrms 1 dcrms 15 s48 .bin 

peouts48.bin 

52m 

peoutzbrmsldcrmsl5s52.bin 

peouts52.bin 

56m 

peoutzbrms 1 dcrms 15s56.bin 

peouts56.bin 

60m 

peoutzbrms 1 dcrms 15s60.bin 

peouts60.bin 

64m 

peoutzbrms 1 dcrms 15s64.bin 

peouts64.bin 

68m 

peoutzbrms 1 dcrms 15s68 .bin 

peouts68.bin 

72m 

peoutzbrms 1 dcrms 15s72.bin 

peouts72.bin 

76m 

peoutzbrms 1 dcrms 15s76.bin 

peouts76.bin 

80m 

peoutzbrmsldcrmsl5s80.bin 

peouts80.bin 


Table 4.2 -MMPEREV/MMPE Output Files 


We can, from a single CW run of MMPE and MMPEREV, make qualitative 
comparisons of the effects of the interface and volume perturbations. 
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Depth (m) Depth (m) 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 


Range (km) 

Figure 4.2 - Transmission Loss (dB re 1 m) with No Perturbation, 
Source: 48 m, Frequency: 250 Hz 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 


Range (km) 

Figure 4.3 - Transmission Loss (dB re 1 m) with Perturbation, 
Source: 48 m, Frequency: 250 Hz 
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When the perturbation was introduced, the structure in Figure 4.2 was disrupted as shown 
in Figure 4.3. The dispersed penetration into the volume is also apparent. In the next few 
sections, we will discuss the post-processing implemented to compute reverberation, 
vertical correlation, and spectral analysis. 

B. POST-PROCESSING 1 - REVERBERATION (MATLAB 
IMPLEMENTATION) 

Recall that in Chapter II we derived the CW reverberant field for a scattering 
patch due to the interface and volume in Eqs. (2.33) and (2.61), respectively. In order to 
arrive at the reverberant field at the receiver due to a co-located or displaced source, we 
had to find the interactions of these fields over the entire scattering area. We will also 
perform broadband and CW computations by utilizing the same output from MMPEREV 
set to run on the various geometrical configurations for the VLA described in the 
previous sub-section. Again, we were able to separate the bottom interface and bottom 
volume analyses because of their different perturbations. The separate analysis also 
afforded closer discrimination of the structure of one independently of the other. 

The source was located at 48 m depth, which was in the vicinity of the center of 
the water column. Each of the array elements (including the 48 m array element) were 
then treated as independent receivers. Each of the output files from MMPEREV and 
MMPE could be paired with the output file for the 48 m source depth binary files for 
computation. Naturally, the 48 m source and 48 m receiver files were the same binary 
file, meaning a monostatic geometry while the other transmitter/receiver pairs reflected a 
vertical bistatic geometry. 

Both perturbed and unperturbed data generated using MMPEREV and MMPE, 
respectively, were analyzed to provide useful comparisons and further analyses. We will 
also show the results for each analysis performed. 

1. Two-way Mono/Bistatic CW Reverberation Analysis 

For the interface, we have from Eq. (2.33) and subsequently in Eq. (2.43), 

55 



(4.13) 


(|P-f ) ' 

This is simply proportional to the magnitude-squared of the forward and backward 
propagated field functions. As the perturbation spectrum and other field constants are 
generated by MMPEREV, the MATLAB post-processing that computes the CW 
reverberation is just the multiplication of the magnitude-squared field functions with 
cylindrical spreading. 


Wn (r) 

Vr 1 4r 

and the interface reverberation loss (RLb) is then 


(4.14) 


RL b (r) = A- 20log 

r 

where A is a constant which accounts for the other parameters in Eq. (2.50). 


(4.15) 


The volume was also treated as described in Eq. (2.66b) and the volume 
reverberation loss may be written as 


RLy(r) = B- 20 log 


- \n{r,z)\yr T (.r,z)\y R (,r,z)\dz 

^ z>z b 


(4.16) 


where B is another constant which accounts for the other terms in Eq. (2.66b). Using now 
the same 16 sets of results from a single MMPEREV and another 16 sets of results from 
the MMPE runs, we post-processed the field function information to arrive at the CW 
interface and volume analyses. We have selected the minimum, center and maximum 
frequencies for analysis. The MATLAB files written to post-process the 
MMPEREV/MMPE output files to arrive at the CW analysis are given in Table 4.3. The 
output from each run of peout2way_bistat_r.m was saved as a MATLAB data file so that 
the data set could be loaded into the MATLAB workspace without having to re-process 
the MMPEREV/MMPE output binary files. The filenames saved are listed in Table 4.5. 
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MATLAB Post-Processing Files 
Description 

Perturbed Data 
Processing Files 

Unperturbed Data 
Processing Files 

Binary file data extractor: Opens both source and 
receiver binary files (output from MMPE/ 
MMPEREV) and all modeling data. Stores 
them on the MATLAB workspace. 

peoutl_bistat_r.m 

peout l_bistat_np_r.m 

Interface and Volume Reverberation loss 
Computation: Extracts field data source and 
receiver files at selected frequencies. 

Computes the two-way reverberation loss for 
the interface and volume as described above. 
Extracts field data at/ mfn ,/ re ™„ and /™*. 

peout2way_bistat_r.m 

peout2way_bistat„np_r.m 


Table 4.4 - MATLAB Filenames Created for CW Reverberation Analysis 


Receiver 

Depth 

Interface 
Perturbed Data 
Ouput 

Volume 
Perturbed Data 
Ouput 

Interface 

Unperturbed Data 
Ouput 

Volume 

Unperturbed Data 
Ouput 

20m 

bistat4820int_r.mat 


bistat4820int_np_r.mat 

bistat4820vol_np_r.mat 

24m 

bistat4824int_r.mat 

bistat4824voLr.mat 

bistat4824int_np_r.mat 

bistat4824vol_np_r.mat 

28m 

bistat4828int_r.mat 

bistat4828vol_r.mat 

bistat4828int_np_r.mat 

bistat4828vol_np_r.mat 

32m 

bistat4832int_r.mat 


bistat4832int_np_r.mat 

bistat4832vol_np_r.mat 

36m 

bistat4836int_r.mat 


bistat4836int_np_r.mat 

bi stat4836vol_np_r. mat 

40m 

bistat4840int_r.mat 

bistat4840vol_r.mat 

bistat4840int_np_r.mat 

bistat4840voLnp_r.mat 

44m 

bistat4844int_r.mat 

■JUMMUUM— 

bistat4844int_np_r.mat 

bistat4844vol_np_r.mat 

48m 

bistat4848int_r.mat 


bistat4848int_np_r.mat 

bistat4848vol_np_r.mat 

52m 

bistat4852int_j\mat 

bistat4852vol_r.mat 

bistat4852int_np_r.mat 

bistat4852vol_np_r.mat 

56m 

bistat4856int_r.mat 


bistat4856int_np_r.mat 

bistat4856vol_np_r.mat 

60m 

unamiiuH 


bistat4860int_np_r.mat 

bistat4860vol_np_r.mat 

64m 



bistat4864int_np„r.mat 

bistat4864vol_np_r.mat 


bistat4868int_r.mat 


bistat4868int_np_r.mat 

bistat4868vol_np_r.mat 

wtEzm 

bistat4872int_r.mat 


bistat4872int_np_r. mat 

bistat4872voI_np_r.mat 

76m 

bistat4876int_r.mat 

^ i j ^ Mi m k I'M 

bistat4876int_np_r.mat 

bistat4876vol_np_r.mat 

80m 

bistat4880int_r.mat 


bistat4880int_np_r.mat 

bistat4880vol_np_r.mat 


Table 4.5 - Processed Output Data Files in CW Analysis 


2. CW Analysis Results 

A typical CW reverberation loss plot shows a 20 log r (two-way cylindrical 
spreading) drop in range. Figures 4.4 and 4.5 show the interface and volume 
reverberation loss for a source depth of 48 m, receiver depth of 40 m and frequency of 
250 Hz. In all cases, the constant scaling parameters A and B were set equal to zero. With 
the same source and receiver depths and at the same center frequency, we were able to 
observe finer structures in the interface reverberation (Figure 4.4) compared to the 
volume reverberation plot (Figure 4.5). This was because multiple incoherent 
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multiplications were performed at all depths and then summed vertically m order to 
arrive at the mean squared RPL at the receiver due to the volume. Note that all data sets 
show very little energy interacting with the interface or volume at very short ranges, i.e. 
< 0.1km. This is due to the wide-angle limitation of the parabolic approximation. 

Source depth:48m, Receiver depth:40m, Frequency: 250Hz 

70 j—-1 I l I l I . i . i ~ i . . 

:::::: — Interface reverberation loss 

in — 20log r drop off 


Range (km) 


Figure 4.4- Interface Reverberation Loss for Two-Way Transmission With Perturbation 


Source depth:48m, Receiver depth:40m, Frequency: 250Hz 


Interface reverberation loss 
20log r drop off _ 


Figure 4.5 - Volume Reverberation Loss for Two-Way Transmission With Perturbation 






















A comparison between perturbed and unperturbed data was also done. Figures 4.6 
and 4.7 show the difference (perturbed RL minus unperturbed RL) to see the fine-scale 
structures due to the introduction of the perturbation. Apparent in these plots are 
repeating structures at about every 0.5 km for ranges beyond about 2 km. 

(Perturbed - Unperturbed) Source Depth: 48m, Receiver depth: 40m, Frequency: 250Hz 



Figure 4.6 - Interface Reverberation Loss Difference Vs Range 
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(Perturbed - Unperturbed) Source depth: 48m, Receiver depth: 40m, Frequency: 250Hz 



Figure 4.7 - Volume Reverberation Loss Difference Vs Range 


Figures 4.8 and 4.9 show the interface and volume reverberation loss, 
respectively, across depths plotted against range for the center frequency of 250 Hz. It 
was apparent that the interface reverberation showed more structure while volume 
reverberation seemed more ‘blotchy’ and revealed little structure. This smudging of data 
was manifest of the multiple incoherent (magnitude) multiplications computed before 
vertically summing the reverberation in depth. What was apparent however, was the 
‘201og r' drop in range in both plots. 
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3. Two-way Mono/Bistatic Time-Domain Reverberation Analysis 

Time-domain analyses for both the interface and volume reverberation were 
explained in Chapter II, Section C. In essence, multiple frequencies were used to simulate 
a pulsed propagation. We will use Eqs. (2.68b) and (2.68c) to formulate the MATLAB 
implementation for the interface and Eqs. (2.70a) to (2.72) for the volume. MMPE was 
also run based on the same 16 setups for the VLA to provide comparison data for 
unperturbed data. Implementation was completed using MATLAB with the files created 


as follows: 


MATLAB Post-Processing Files 
Description 

Perturbed Data 
Processing Files 

Unperturbed Data 
Processing Files 

Binary file data extractor: Opens both source 
and receiver binary files (output from 
MMPE/MMPEREV). Extracts all modeling 
parameters and stores them on the 

MATLAB workspace. 

Peoutl_bistat.m 

peoutl_bistat_np.m 

Interface and Volume Reverberation 

Computation: Extracts from binary files the 
field function data for source and receiver 
and computes the two-way reverberation 
loss for die interface and volume as 
described above. 

Peout2way_J)istat.m 

peout2 way_b i s tat_np. m 


Table 4.6 - MATLAB Filenames Created for Time-Domain Reverberation Analysis 

The output from each run of peout2way_bistat.m was saved as a MATLAB data file so 
that the data set could be loaded into the MATLAB workspace without having to re¬ 
process the MMPEREV/MMPE output binary files. These are listed in Table 4.7. 

Additionally, peoutl_bistat.m and peout2way_bistat.m were modified into 
function form in order to allow unattended batch processing of the transmitter/receiver 
pairs. MATLAB routines created were: master_run.m, peoutl_bifn.m, 
peout2way_biintfn.m, and peout2way_bivolfn.m, where the interface and volume 
processing in the original peout2way_bistat.m was split up into two files for processing 
efficiency. 
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Receiver 

Depth 

Interface 
Perturbed Data 
Ouput 

Volume 
Perturbed Data 
Ouput 

Interface 
Unperturbed 
Data Ouput 

Volume 
Unperturbed 
Data Ouput 

20m 

bistat4820int.mat 

bistat4820vol.mat 

bistat4820int_np.mat 

bistat4820vol_np.mat 

24m 

bistat4824int.mat 

bistat4824vol.mat 

bistat4824int_np.mat 

bistat4824vol_np.mat 

28m 

bistat4828int.mat 

bistat4828vol.mat 

bistat4828int_np .mat 

bistat4828vol_np.mat 

! 32m 

bistat4832intmat 

bistat4832vol.mat 

bistat4832int_np.mat 

bistat4832vol_np.mat 

36m 

bistat4836int.mat 

bistat4836vol.mat 

bistat4836int_np.mat 

bistat4836vol_np.mat 

40m 

bistat4840int.mat 

bistat4840vol.mat 

bistat4840int_np.mat 

bistat4840vol_np.mat 

44m 

bistat4844int.mat 

bistat4844vol.mat 

bistat4844int_np.mat 

bistat4844vol_np.mat 

48m 

bistat4848int.mat 

bistat4848vol.mat 

bistat4848int_np.mat 

bistat4848vol_np.mat 

52m 

bistat4852int.mat 

bistat4852vol.mat 

bistat4852int_np.mat 

bistat4852vol_np.mat 

56m 

bistat4856int.mat 

bistat4856vol.mat 

bistat4856int_np.mat 

bistat4856vol_np.mat 

60m 

bistat4860int.mat 

bistat4860vol.mat 

bistat4860int_np.mat 

bistat4860vol_np.mat 

64m 

bistat4864int.mat 

bistat4864vol.mat 

bistat4864int_np.mat 

bistat4864vol_np.mat 

68m 

bistat4868int.mat 

bistat4868vol.mat 

bistat4868int_np.mat 

bistat4868vol_np.mat 


bistat4872int.mat 

bistat4872vol.mat 

bistat4872int_np.mat 

bistat4872vol_np.mat 


bistat4876int.mat 

bistat4876vol.mat 

bistat4876int_np.mat 

bistat4876vol_np.mat 

80m 

bistat4880int.mat 

bistat4880vol.mat 

bistat4880int_np.mat 

bistat4880vol_np.mat 


Table 4.7 - Processed Output Data Files in Time-Domain Analysis 


4. Time-Domain Analysis Results 

Having discussed the implementations, we now show some results obtained. First, 
we looked at the general structure of interface and volume reverberation in the time- 
domain at a particular depth chosen arbitrarily as 40 m. Showing very much the same 
characteristics, Figures 4.10 and 4.11 depict the time-domain interface and volume 
reverberation loss, respectively. The figures provide an understanding of the structure of 
the reverberation return of a pulsed transmission with returns lasting for the entire extent 
of the calculations both in the interface and volume. 
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Figure 4.10 - Time-Domain Interface Reverberation Loss 




































Time-domain data sets were also used to provide useful comparisons with 
unperturbed data sets. To do so, the time axis was ‘converted’ to a range axis by use of 
the reference sound speed, c 0 , such that range = c 0 x time . This produced a range axis 

from 0 to 5 km. The data comparison then proceeded along the same path as that in the 
CW analysis in the previous section. We first show sample time-series computed 
reverberation-range plots for the interface and volume. Shown in Figures 4.12 and 4.13 
are the time-domain reverberation loss plots equivalent to Figures 4.4 and 4.5 for the CW 
analyses. Again, the ‘201og r’ drop off was apparent. 


Source depth: 48m, Receiver depth: 40m 



Figure 4.12 -Time-Domain Converted Interface Reverberation Loss 
(dB re 1 m) Vs Range, With Perturbation 
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Source depth: 48m, Receiver depth: 48m 



Figure 4.13 - Time-Domain (Converted to Range) 
Volume Reverberation Loss Vs Range, With Perturbation 


Now looking across depths, we observe multipath structures shown in Figures 
4.14 and 4.15 due to both interface and volume perturbations, as expected. Spatial 
disposition of each source was important in observing the arrival time of the first returns, 
most apparent at the receiver element at 20 m depth, while arrivals were earlier at the 
receiver element at 80 m depth. This was due simply to the closer proximity of the deep 
receiver to the bottom thereby receiving returns earlier than the shallow receivers. No 
other structures were apparent from these two figures. 
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Source Depth: 48m 



Time (sec) 


Figure 4.14 - Interface Reverberation Loss Across Depth Vs Time 


Source Depth: 48m 



Time (sec) 


Figure 4.15 - Volume Reverberation Loss Across Depth Vs Time 
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5. MATLAB Data Analysis 

As a final note on the work done in both data and signal analyses, several 
MATLAB files were also written to ensure that the data analysis was complete. These are 
found in Table 4.8. In order to derive further value, vertical correlations and spectral 
analyses were performed. These are discussed in the next sections. 


MATLAB Post-Processing Files Description 

MATLAB Files 

Range-Series Data Processor: Using data obtained from 
the output file created by peout31_pre_bistat_r.m and 
peout31_pre_bistat_r_np.m respectively, this routine 
computes and displays reverberation loss and 
compares these RLs between perturbed and 
unperturbed range-series data. 

range_series_datanalysis. m 

Time-Series Data Processor: Using data obtained from 
the output file created by peout31_pre_bistat_r.m and 
peout31_pre_bistat_r_np.m respectively, this routine 
computes and displays reverberation loss and 
compares these RLs between perturbed and 
unperturbed range-series data. 

time_series_datanalysis.m 


Table 4.8 - MATLAB Routines Created for Data Analyses and Comparisons 


C. POST-PROCESSING 2 - VERTICAL CORRELATION & PEAK 
CORRELATION IN CW AND TIME-DOMAIN 

Having completed CW and time-domain analyses, correlation between the 
monostatic (source and receiver co-located at 48 m) reverberation data set with 
reverberation data sets from all other receiver depths was performed. The data files saved 
from the CW (Table 4.5) and time-domain (Table 4.7) analyses were used to compute the 
vertical correlations for both CW (range correlation) and time-domain (temporal 
correlation). MATLAB routines were written to perform the vertical correlations. All 
vertical correlations are computed relative to a source depth of 48 m. Thus, the 
autocorrelation would occur at the receiver depth of 48 m, giving the theoretical 
maximum correlation over depth. These were done for both the interface and volume 
with and without perturbations. Selected results are shown in this section. 
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1. Vertical Correlation in Range & Peak Vertical Correlation 

We first looked at CW vertical correlation. Table 4.9 lists the MATLAB files 
created to compute the vertical correlations. We have selected a source at 48 m and at the 
center frequency of 250 Hz for this analysis. With the perturbation included, Figure 4.16 
shows the vertical correlation for the interface while Figure 4.17 shows that in the case of 
the volume. It should be noted that due to the incoherent processing of the two-way 
pressure fields, the first 0.5 km does hot accurately reflect the true reverberation loss. 
Hence, the first 0.5 km has been truncated prior to performing the volume vertical 


correlations shown in Figure 4.17. 


MATLAB Post-Processing Files 
Description 

MATLAB Files 
(Perturbed) 

MATLAB Files 
(Unperturbed) 

Correlation Computation: Extracted data from 
Table 4.5 and computed vertical and peak 
vertical correlations with the data at source 
depth, 48m. 

Peout3_bistat_r.m 

peout3_bistat_r_np .m 

Data Processing: Extracted data from Table 4.5 
and stored workspace variables into an 
output data file for future analyses. 

Peout3 l_pre_bistat_r.m 

peout31 _pre_bistat_r_np.m 

Correlation Computation: Computes vertical 
and peak vertical correlations as in 
peout3_bistat_r.m and 
peout3_bistat_r_np.m respectively, but uses 
output file created by 
peout31_pre_bistat_r.m and 
peout3 l_pre_bistat_r_np.m respectively. 

Peout3 l_bistat jr.m 

peout3 l_bistat_r_np.m 

Output filename created by 
peout3 l_pre__bistat_r.m and 
peout3 l_pre_bistat_r_np.m respectively 

Bistat_depmg_data.mat 

bistat_depmg_np_data.mat 


Table 4.9 - MATLAB Filenames Created for Computing Vertical Correlations 
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Relative Depth (m) Relative Depth (m) 


Source Depth: 48m, Frequency: 250Hz 



Range Lag (km) 

Figure 4.16 - Vertical Correlation of Interface Reverberation Loss 
in Range Vs Relative Depth 


Source Depth: 48m, Frequency: 250Hz 



Range Lag (km) 


Figure 4.17 - Vertical Correlation of Volume Reverberation Loss 
in Range Vs Relative Depth (first 0.5 km truncated) 
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It was apparent that finer scale structures were present in the interface, but the volume 
showed strong correlation throughout depth for much of the computations in range up to 
about 3 km due to the smoothing effect of the incoherent processing of the CW signal. 

Peak correlation values were also extracted from the above vertical correlation 
structures in order to see how rapidly the signal decorrelated over depth. As would be 
expected from the above analysis, the volume would tend to stay more correlated over 
depth while the interface showed more decorrelation over depth. However, even the 
interface reverberation loss remained quite correlated with minimum correlation values 
down only about 1 dB. This is also due presumably to the incoherent processing involved 
in the derivation of the mean squared RPL. Figures 4.18 and 4.19 show a comparison of 
the peak vertical correlations between perturbed and unperturbed data. While it is 
apparent that the volume reverberation of the CW signal tends to decorrelate slower than 
the interface in both scenarios, the addition of the perturbations does not appear to affect 
the vertical correlations significantly. 


Source Depth: 48m, Frequency: 250Hz 

-30| i i i i i i i r 



|_ I _l_ I _i_i_1_ I _i_I 

-1.6 -1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 

Peak Correlation (dB) 

Figure 4.18 - Peak Vertical Correlations (dB) of 
Perturbed Data Over Relative Depth 
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Source Depth: 48m, Frequency: 250Hz 
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-1.6 -1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 

Peak Correlation (dB) 

Figure 4.19 - Peak Vertical Correlations (dB) of 
Unperturbed Data Over Relative Depth 

2. Vertical Temporal Correlation & Peak Vertical Temporal 
Correlation 

As with the CW signals, MATLAB files were created to compute the vertical 
correlations in the time-domain for the broadband signals. These are listed in Table 4.9. 
From the results of the vertical temporal correlations and peak vertical temporal 
correlations, we obtain Figures 4.20 through 4.23. 

Apparent from Figures 4.20 and 4.21 are finer scale structures seen in both the 
interface and volume reverberation level. Of importance was the fact that the coherent 
convolution done in arriving at the time-domain analysis preserved the coherent time 
structures, in contrast to the CW analyses. 
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MATLAB Post-Processing Files 
Description 

MATLAB Files 
(Perturbed) 

MATLAB Files 
(Unperturbed) 

Correlation Computation: Extracted data from 
Table 4.7 and computed vertical and peak 
vertical correlations with the data at source 
depth, 48m. 

peout3_bistat.m 

peout3_bistat_np.m 

Data Processing: Extracted data from Table 4.7 
and stored workspace variables into an 
output data file for future analyses. 

peout3 l_pre_bistat.m 

peout3 l_pre_bistat_np.m 

Correlation Computation: Computes vertical 
and peak vertical correlations as in 
peout3__bistat.m and peout3_bistat_np.m 
respectively, but uses output file created by 
peout31_pre_bistat.m and 
peout3 l_prej)istat_np.m respectively. 

peout31_bistat.m 

peout3 l_bistat_np.m 

Output filename created by 
peout31_pre_bistat.m and 
peout3 l_pre_bistat_np.m respectively 

bistat_deptime_data.mat 

! 

bistat_deptime_np _data.mat 


Table 4.10 - MATLAB Filenames Created for Computing Vertical Correlations 


Source Depth: 48m 



Time Lag (sec) 



Figure 4.20 - Vertical Temporal Correlation of 
Interface Reverberation Loss Vs Relative Depth 
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Source Depth: 48m 



Time Lag (sec) 

Figure 4.21 - Vertical Temporal Correlation of 
Volume Reverberation Loss Vs Relative Depth 

Peak vertical temporal correlations were extracted for the interface and volume 
relative to a 48m source, and the results are displayed in Figure 4.22. Similar analysis 
was performed on data computed without the perturbations and these results are 
presented in Figure 4.23. With the same horizontal scales, Figures 4.22 and 4.23 could be 
compared to see the effects of introducing perturbations. The perturbations appear to 
affect the volume reverberation loss vertical correlation more than the interface 
reverberation loss, although the effect is minor in both cases. For both perturbed and 
unperturbed environments, the vertical 3dB decorrelation length scale is on the order of a 
few meters. 
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Source Depth: 48m 



Figure 4.22 - Peak Vertical Temporal Correlations 
of Perturbed Data 


Source Depth: 48m 



Figure 4.23 - Peak Vertical Temporal Correlations 
of Unperturbed Data 


75 










































D. POST-PROCESSING 3 - SPECTRAL ANALYSES 


In this final section, we investigate the spectral content of the reverberation 
pressure obtained for both the CW and broadband signals. There was a need to convert 
the time-domain into range using the reference sound speed, Co, like the conversion done 
in Section 3.B in the time-domain analysis. As before, comparison with unperturbed data 
was also done. MATLAB routines were again written and tabulated as follows. 


MATLAB Post-Processing Files 
Description 

MATLAB Files 

Computes the Fourier Transform of the CW 
data obtained from the output file created by 
peout3 l_pre_bistat.m and 
peout31_pre_bistat_np.m respectively. This 
analysis has built within itself a comparison 
with unperturbed data already. 

range_series_signalysis.m 

Computes the Fourier Transform of the 

Converted Time-Domain data obtained from 
the output file created by 
peout31_pre_bistat.m and 
peout31_pre__bistat_np.m respectively. This 
analysis has built within itself, a comparison 
with unperturbed data already. 

time_series_signalysis.m 


Table 4.11 - MATLAB Routines Created for Spectral Decomposition 


Assuming a monostatic geometry, we recall from the analyses in Chapter II and 
from Eqs. (4.15) and (4.16) that the interface reverberation loss may be defined as 

RI^ = 2 TL b + constants , (4-17) 

where 


TL b =101og|p_(r)| 2 =201og|^(r)|-101ogr 
since, by definition, 

\p(r)\ = -j=\y(r)\. 

Hence 

RLh °c 20 log|^ (r)| 2 - 20 log r. 

We define the range-reduced reverberation loss as 


(4.18) 


(4.19) 


(4.20) 
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(4.21) 


RL bjeduced =RL b + 20log r oc 20 \og\y/ b (r)| 2 . 

Now we assume that we can represent the field function in terms of a spectrum, 
i.e. 

\n<.ri~\w(K)e‘ K 'dK. (4 ' 22) 

It is possible then that the spectral content of the field can be related to the statistical 
characteristics of the perturbations and can thus be extracted with the use of some signal 
processing tools. In other words, the spectrum W(K) may be related to the perturbation 
spectra. 

The main approach was the use of discrete Fourier transforms (DFT) in an 
attempt to extract the spectral components of the reverberation data. They are: 

• Signal Analysis 1: Power Spectral Density - DFT of the magnitude squared 
range-reduced reverberation (i.e. r 2 x|p_| 2 ). 


• Signal Analysis 2: Power Ratio Spectral Density - DFT of the magnitude 
squared of the ratio of the reverberation (perturbed divided by unperturbed), 


i.e., 


\P~ perturbed 


P-unperturbed 


2 * 


Recall that the typical reverberation loss observed from the range series (CW) 
data for the interface and volume were presented earlier in this chapter in Figures 4.4 and 

4.5 respectively. The time-domain reverberation loss data was correspondingly shown in 
Figures 4.10 and 4.11. The difference in RL (ratio) between the perturbed and 
unperturbed data (CW analysis) for the interface and volume were also shown in Figures 

4.6 and 4.7. Figures 4.24 and 4.25 show the CW reverberation loss while Figures 4.26 
and 4.27 show the time-domain (converted to range) reverberation loss. 
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Source depth: 48, Receiver depth: 40m, Frequency: 250Hz 



Figure 4.24 - Interface Reverberation Loss vs 
Range Using Range Series Data 


Source depth: 48, Receiver depth: 40m, Frequency: 250Hz 



Figure 4.25 - Volume Reverberation Loss vs 
Range Using Range Series Data 

When using time-series data to perform this analysis, time was converted to range 
in order to relate length scales of the reverberation with perturbation scales. It is perhaps 
interesting to note that the time-series data showed no appreciable differences between 
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the perturbed and unperturbed data sets, in contrast to the CW range-series data displayed 
above. The reason for this is currently unknown. 


(Time Series) Source depth: 48m, Receiver depth: 40m 



Figure 4.26 - Interface Reverberation Loss (dB re 1 m) vs 
Range Using Time Series Data 

(Time Series) Source depth: 48m, Receiver depth: 40m 
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1. Signal Analysis 1 - Power Spectral Density 

The magnitude-squared of the range-reduced perturbed reverberation was 
analyzed using the DFT, giving 

PSD = DFT{\p_\ 2 x r 2 } (4.23) 

where p- is the reverberation field of the interface or volume. Figures 4.28 and 4.29 show 
the normalized power spectrum of the range-reduced interface and volume 
reverberations, respectively. Note that these plots are presented with the wavenumber 
axis displayed on a logarithm scale in order to elucidate exponential behavior. While 
structurally the same, i.e., lower wavenumbers predominating the spectra, strong fine 
scale structures were evident in both plots for the perturbed data. This alludes to the 
presence of perturbations which may be extracted from the spectra. In the volume plot, 
Figure 4.29, the perturbed data distinctly showed less drop-off beyond a wavenumber of 
roughly 0.2 m' 1 corresponding to wavelengths smaller than approximately 30 m. 


Source depth: 48m, Receiver depth: 40m, Frequency: 250Hz 



Wavenumber, (nrf 1 ) 


Figure 4.28 - Normalized Power Spectrum of 
Range-Reduced Interface Reverberation (Range-Series) 
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Source depth: 48m, Receiver depth: 40m, Frequency: 250Hz 



Wavenumber, (m -1 ) 

Figure 4.29 - Normalized Power Spectrum of 
Range-Reduced Volume Reverberation (Range-Series) 

Figures 4.30 and 4.31 show the results for the broadband signal (after converting 
time to range) for the interface and volume. Both fine and large scale structures are 
observed. While repeating structures were mostly preserved when the perturbation was 
included, the longer scale perturbations may have augmented the spectra. 
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Signal Analysis 2 - Power Ratio Spectral Density 


The quotient of the perturbed over the unperturbed reverberation data was 
magnitude-squared to provide a comparison of the levels. This was then Fourier 
transformed in order to reveal any special spectral content. The power ratio spectral 
density (. PRSD ) was then defined as 


PRSD = DFT 



P-perturbed 

2 ' 


P -unperturbed 

T * 


(4.24) 


Figures 4.32 and 4.33 show the spectral components of the ratio between the perturbed 
and unperturbed data sets for the CW signal. While it was difficult to make any apparent 
observations, the interface reverberation data showed a rather flat spectrum while the 
volume reverberation data showed a roll-off at higher wavenumbers. The physical 
justification for this is, as yet, unclear. 


(IPerturbedp/jUnperturbedl 2 ) Source depth: 48m, Receiver depth: 40m, Frequency: 250Hz 














































(|Perturbed| 2 /|Unperturbed| 2 ) Source depth: 48m, Receiver depth: 40m, Frequency: 250Hz 
Oi-1-1—*—; rrm -:- t —:— i rrm -:-r—r—" 



Wavenumber, (m _1 ) 


Figure 4.33 - Normalized Spectrum of the Volume Power Ratio 
Using Range-Series Data 

Finally, we computed the power ratio spectral density based on the time-series 
data, as shown in Figures 4.34 and 4.35 for the interface and volume, respectively. These 
plots were consistent with the above CW analysis for the ratio of the perturbed to the 
unperturbed data. The interface showed an almost flat spectra while the volume showed a 
drop-off with increaasing wavenumber. 
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(IPerturbed^/IUnperturbedl 2 ) Source depth: 48m, Receiver depth: 40m 



Figure 4.35 - Normalized Spectrum of Volume Power Ratio 
Using Time-Series Data 
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V. CONCLUSIONS & RECOMMENDATIONS 


A theoretical model and corresponding numerical simulation for shallow water 
reverberation has been developed for spatial and sound speed perturbations of the bottom 
interface and volume, in the context of the parabolic approximation. The main focus in 
this thesis has been the use of the MMPE model for predicting the forward and back 
propagations, and then developing post-processing routines in MATLAB to compute the 
reverberation loss and to perform signal analyses on the derived reverberation. 

The development of statistical models to introduce interface and volume 
perturbations into the MMPE model were also completed, noting that density variations 
to the volume have yet to be included to arrive at a more comprehensive model. This is, 
however, not envisaged to be difficult and will be left for follow-on work. The models 
presented here were tested and found to be effective in modeling the perturbations typical 
of different types of bottoms and sediment compositions. An rms value of 1 m interface 
perturbation and a 15 m/s sound speed variation were used. Most apparent was the effect 
these perturbations had on a one-way propagation of the MMPE model. It showed 
diffused bottom penetration extending much farther in range with less noticeable effects 
in the water column. 

Various analyses were conducted for both perturbed and unperturbed 
environments in order to establish baselines for which to see the effects of perturbation. 
This was done in computing the reverberation losses across a VLA based on a single 
source located at range, r = 0m, and at 48 m depth. Arrival structures were clearly 
discernible due to travel time differences, the shallowest receivers showing later arrival 
times. Reverberation levels consistently showed ‘201og r' drop-offs regardless of the 
presence of perturbation. Correlation across depth was performed and showed that the 
interface decorrelated more rapidly than for the volume. This was also true of the 
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scenario without perturbations. Thus, the perturbations were not seen to affect correlation 
significantly. 

The power spectrum of the reverberation and of the ratio of the perturbed to the 
unperturbed data were computed. In the power spectral analysis of the volume 
reverberation, showed less drop-off for the perturbed data. However, both interface and 
volume-perturbed plots showed small and large-scale deviations from the unperturbed 
plots, providing evidence of the perturbations. Broadband analysis showed that repeating 
structures were preserved while shorter scale perturbations may have augmented the 
spectra. The PRSD showed a rather flat spectrum for the interface reverberation data 
while volume reverberation showed roll-off at higher wavenumbers. This was consistent 
in both the CW and broadband analyses but the justification is not presently known. 

In essence, we have analyzed the influence of interface and volume perturbations 
on CW and broadband reverberation structure in shallow water. The impact of the 
interface perturbations cannot be really concluded at this time since the spectral analyses 
did not provide any observable characteristics. The volume perturbation spectrum 
appeared more amenable to extraction via signal processing and should be explored 
further. 

With the conclusion of this thesis, recommendations for future work are: 

• to generalize the Fortran implementation of the perturbation models with truly 
random perturbations by incorporating a ‘seed’; 

• to develop and implement density fluctuations into the sediment volume and 
reverberation theory; 

• and to perform second and third order signal analyses on the reverberation 
data in an attempt to extract information on the perturbation spectra. 
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